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A winter  survey  was  conducted  in  Waconda  Bay,  Lake  Chickamauga.  The 
Bay  is  the  receiving  water  for  TNT-manufacturing  wastes  from  Volunteer 
Army  Ammunition  Plant  located  northwest  of  Chattanooga,  Tennessee. 


Field  work  was  carried  out  in  two  phases.  The  first,  completed  in 
'December,  determined  flow  and  dispersion  patterns  in  Waconda  Bay,  munitions 
levels  in  the  lake  water  and  sediments  and  other  limnological  parameters 
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20.  Abstract  (cont. ) 

/ ^of  the  system.  The  Phase  I surv^  detennined  that  munitions  levels  were 

sufficiently  high  to  justify  an  intensive  biological  study  at  a later 

date  and  provided  data  for  locating  sampling  sites  covering  the  full 
range  of  concentrations. 


The  second  part  of  the  investigation  was  completed  during  March  and 
examined  tHe  effects /of  TNT  wastes  from  VAAP  on  the  periphyton  and  macro- 
invertebrate communities,  A rigqrous  supplementary  study  determined  the 
statistical  variability  ociiur’ring  in  periphyton  colonizing  glass  slides. 

Utilizing  the  population  density  and  community  structure  of  these 
selected  biological  compartments,  estimates  were  made  of  effect  levels 
from  munitions  waste.  /In  the  natural  substrate  no  effects  could  be 
discerned  in  either  number  of  species,  population  density, or  community 
structure  for  both  periphyton  or  macroinvertebrates  at  locations  3,000 
feet  downbay  from  tbe  point  of  waste  entry  where  total  TNT  concentration 
averaged  less  than  2'5'ug/l.  During  the  30  day  exposure  period  for 
artificial  substrates,  samples  of  VAAP  effluent  were  collected  and 
• ’analyzed  for  TNT,  Concentration  of  5 munitions  compounds  averaged  575  ug/1 
and  ranged  as  high  as  2,370  wg/1.  In  the  0-3,000  foot  zone,  definite 
•effects  were  noted  in  the  periphyton  and  macroinvertebrates  colonizing 
artificial  substrates.  Effects  were  manifested  by  shifts  in  population, 
number  of  species,  and  community  structure. 


the  results  of  the  study  indicate  no  environmental  effect  on  peri- 
phyton or  macroinvertebrates  from  a complex  TNT  manufacturing  effluent 
at  a total  munitions  concentration  of  less  than  25  ng/1.  Definite 
effects  were  noted  during  a period  wheo^total  munitions  averaged  500- 
60(V^gf^with  a peak  of  over  2,000(^pg^  Concentrations  believed  to  be 
in  the>ange  of  50-100  ug/1  produced  minimal  effects.. 
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EXECUTIVE  SUMMARY 


Field  surveys  were  conducted  during  December,  1976,  and  March,  1977, 
in  Waconda  Bay  to  deteniiine  the  distribution  and  effect  on  periphyton  and 
macroinvertebrates  of  TNT  effluent  discharged  from  Volunteer  Army  Aiiiiiuni tion 
Plant.  The  study  was  an  extension  of  a 1975  sunnier  effort  to  detennine  no- 
effects levels  for  munitions  compounds  in  fresh  water. 

The  study  was  divided  into  two  phases.  Phase  I,  completed  in  December, 
defined  the  movement  of  the  effluent  in  the  bay  to  aid  in  selection  of  chem- 
ical and  biologic  sampling  sites  for  the  more  extensive  Phase  II  work.  A 
dye  tracer  study  utilizing  Rhodamine  B,  in  situ,  and  total  munitions  measure- 
ments were  utilized  in  Phase  1.  Results  showed  that  dispersion  of  effluent 
into  the  Bay  was  most  influenced  by  reservoir  stage  changes.  A rising  stage 
had  the  effect  of  holding  the  effluent  in  the  upper  part  of  the  Bay  while  a 
drawdown  period  moved  effluent  rapidly  downbay. 

During  a one-week  intensive  sampling  period,  total  munitions  concen- 
tration was  highest  in  the  first  3,000  feet  from  the  waste  outfall.  Mean 
total  levels  were  as  high  as  397  pg/l.  The  principal  munitions  compounds 
near  the  discharge  were  2,4-DNT,  2,6-DNT,  and  2,4,6-TNT.  Analysis  of  data 
from  this  first  phase  shows  that  pollutant  concentration  in  the  initial  3,000 
feet  was  about  0.5  to  1 times  the  effluent  concentration.  Levels  diminished 
rapidly  the  next  1,600  feet  to  0.01  to  0.05  times  that  of  the  effluent.  Based 
on  results  of  the  December  study,  stations  for  Phase  11  were  grouped  in  high, 
medium,  and  low  concentration  zones.  Chemical  and  biological  characterizations 
were  deteniiined  for  water  and  sediments.  Low  stage  conditions  which  continued 
to  prevail  during  March  (Phase  II)  caused  a channel  to  form,  extending  about 
3,000  feet  beyond  the  outfall. 

Daily  munitions  analysis  during  the  first  week  in  March  showed  that 
highest  levels  occurred  consistently  in  the  channel.  Average  total  munitions 
concentration  in  this  area  was  78  ng/l  with  the  principal  components  being 
2,4-DNT,  2,6-DNT,  and  2,4,6-TNT.  Total  munitions  concentration  in  the  effluent 
during  this  time  based  on  average  24  hour  composites,  was  89  ug/l.  It  is 
apparent  that  little  dilution  was  occurring  in  the  channel.  Concentrations  of 
munitions  compounds  dropped  significantly  beyond  the  channel.  Relative  to 
other  compounds  only  three  of  the  individual  Bay  samples  showed  2,4-DNT  and 
five  samples  contained  1,3-DNR.  Only  1,3,5-TNB  showed  increased  levels 
downbay.  It  appears  that  this  compound  is  a breakdown  product  of  2,4,6-TNT. 

No  1,3,5-TNB  was  detected  in  the  plant  effluent. 

In  early  March  total  nitrogen  in  the  channel  averaged  14.2  mg/1  with 
approximately  79  percent  occurring  as  NO3-N.  Concentrations  downbay  showed  a 
gradual  reduction.  The  N levels  in  the  medium  zone  averaged  5.6  mg/1  with 
nitrate  occurring  again  in  the  same  proportion.  Outer  bay  stations  reflected 
a mean  nitrogen  content  of  4 mg/1  with  NO3-N  concentration  as  the  most 
significant  form.  Other  constituents  such  as  chloride,  total  hardness,  and 
sulfate  followed  a s imi lar pattern  of  decreasing  concentrations  downbay. 
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Foi'  a period  of  one  month  while  artificial  substrates  were  exposed 
in  Waconda  Bay,  24-hour  composite  samples  were  collected  of  VAAP  effluent. 
Results  show  that  average  levels  duriiu]  the  incubation  of  bioloyical  samplers 
were  higher  than  during  the  early  March  field  survey.  Total  munitions  in  the 
channel  pj'obably  were  nea»'  2,000  ug/1  at  a maximum  and  averaged  500  ug/l. 
Organisms  in  the  intermediate  zone  just  beyond  the  channel  were  likely  exposed 
to  a concentration  of  50  - 100  ug/l. 

Periphyton  and  macroinvertebrates  were  the  biological  components 
selected  for  stutty.  Of  these,  periphyton  analysis  provided  the  most  useful 
information.  Community  structure  was  examined  following  a .10-day  colonization 
of  standard  glass  microscope  slides.  Suspension  units  containing  up  to  8 
peri  phytometers  weiv  utilized  during  incubation.  A correlative  study  examined 
between  slide  and  within  slide  variability  using  population  size  and  species 
diversity  as  variables. 

The  periphyton  population  density  was  highly  correlated  with 
concentrations  of  VAAP  effluent.  Population  sizes  expressed  as  no/miir  were 
all  loss  than  1,000  in  the  narrow  channel.  As  the  channel  widened  where  down- 
bay  water  mixed  and  diluted  the  waste,  populations  increased  two  orders  of 
magnitude  and  remained  at  this  level  throughout  Waconda  Bay.  Treatment  of  the 
data  by  the  Pinkham-Pearson  Similarity  Index  shows  this  relationship  among 
stations.  Levels  of  chlorophyll  a and  organic  biomass  add  support  to  this 
latter  conclusion.  Chlorophyll  a on  colonized  slides  at  the  various  stations 
show  a distribution  consistent  with  population  counts.  Organic  biomass 
follows  a simila)’  pattern.  However,  Shannon-Weavet'  diversity  did  not  coi'relate 
with  munitions  concent)'ations . This  index  appears  to  be  of  limited  value  in 
situations  where  the  main  environmental  response  is  a change  in  population 
liens  i ty. 

The  statistical  analysis  indicates  that  a principal  source  of  variability 
in  estimating  either  diversity  or  population  size  of  periphyton  colonizing 
artificial  substrates  is  the  difference  among  direct  microscopic  field  counts 
within  individual  slides.  Unegual  distribution  of  organisms  which  occurred 
during  the  drying  of  the  aliguot  on  the  glass  coverslip  produced  this  vari- 
ability. Locations  near  the  edge  of  the  coverslip  had  lower  populations  than 
those  near  the  center.  Tl\e  tvsult  of  an  uneven  distribution  of  diatoms  is  a 
masking  of  aliguot  and  slide  to  slide  variation.  Tt\e  eri’or  due  to  drying  on 
the  coverslip  must  be  reduced  before  a full  understanding  of  periphyton  vari- 
ab i 1 i ty  can  be  a tta i ned . 

Macwinvertebrate  distribution  in  sediments  did  not  reveal  clearly 
defined  relationships  to  TNT  waste  products.  Low  populations  encountered 
because  of  the  cold  season  ai'e  consideivd  as  the  principal  reason.  Midges 
which  are  an  abundant  fonii  in  Waconda  Bay  colonize  mainly  during  the  suniner 
but  are  sparse  during  the  winter.  Species  dominating  in  the  March  survey 
were  iVpcladius  sp.,  I'pejojanypus^  spp.  and  a f.hmnuwuis  sp. 
dominateif  ilftKe  channeT  nearest  to  the  outfan  anJmay  have  represented  a 
form  tolerant  to  the  waste. 
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Chironotnids  also  dominated  Hester-Dendy  artificial  substrate  samplers. 
The  data  show  a more  defined  relationship  among  stations  than  is  apparent 
with  sediments.  Even  so  the  definition  is  not  clear  cut.  There  appears  to 
be  maximum  impact  in  the  3,000  foot  channel  section  with  a suggestion  of 
transition  zones  in  the  outer  bay  areas.  Overall  the  results  suggest  that 
macroinvertebrates  are  not  as  useful  a parameter  to  study  TNT  effects  during 
a winter  sampling  effort  as  are  periphyton. 

It  should  be  noted  that  the  observed  correlations  between  munitions 
concentrations  and  biologic  responses  do  not  establish  precise  cause  and 
effect  relationships.  The  complex  effluent  at  VAAP  contains  many  compounds 
that  vary  in  concentration  with  regard  to  the  specific  munitions-unique 
compounds  measured.  These  compounds  may  be  factors  also  in  the  effects 
observed. 

Nevertheless,  based  on  the  winter  survey,  no  biological  effects  were 
discerned  when  total  munitions  concentration  was  less  than  25  ug/l.  Definite 
effects  were  noted  when  munitions  levels  averaged  500  - 600  yg/1  with  a peak 
of  around  2,000  yg/1.  Minimal  effects  were  observed  at  munitions  levels  in 
the  range  of  50  - 100  pg/1.  These  data  suggest  a no-effects  level  of  25  yg/1 
for  a complex  TNT  effluent. 
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INTRODUCTION 


The  Volunteer  Army  Ammunition  Plant  (VAAP)  is  situated  on  7,300  acres 
northwest  of  Chattanooga,  Tennessee  and  manufactures  trinitrotoluene  (TNT) 
on  a contractor-ope'-ated  basis.  Wastewater  from  VAAP  drains  northward  into 
a series  of  treatment  lagoons  and  is  discharged  into  the  head  of  Waconda 
Bay  after  undergoing  pH  adjustment  with  lime.  Figure  1 shows  Waconda  Bay 
in  relationship  to  Harrison  Bay  and  Lake  Chickamauga. 

The  purpose  of  this  study  was  to  determine  the  effects  of  munitions 
wastes  on  periphyton  and  macroinvertebrates  during  the  cold  season.  The 
investigation  was  dual-phased  initiated  by  a dye  study  coupled  with  abbreviated 
environmental  sampling.  The  main  objective  of  the  Phase  I study  was  to 
determine  if  sufficient  concentrations  of  munitions  wastes  were  present  in 
Waconda  Bay  to  justify  more  intenstive  biologic  work  and,  if  so,  to  select 
sampling  sites  so  that  the  full  range  of  munitions  concentrations  could  be 
evaluated. 

The  Phase  II  work  included  more  intensive  biologic  sampling  of  peri- 
phyton and  macroinvertebrate  communities  utilizing  both  artificial  and  natural 
substrates.  This  work  was  originally  scheduled  for  early  February  but  was 
postponed  to  early  March  because  Waconda  Bay  was  frozen  in  January  and  February. 
A supplementary  study  was  included  in  Phase  II  to  determine  the  statistical 
variability  occurring  in  population  density  and  community  structure  of  peri- 
phyton colonizing  glass  slides  as  a function  of  sampling  and  processing 
technique. 

Wastewater  and  sediment  characterization  was  conducted  throughout  both 
segments  of  the  field  study.  Plant  effluent  samples  were  also  collected  and 
analyzed  for  TNT  and  associated  transformation  products  for  a one-month  period 
following  the  second  phase  field  survey.  This  coincided  with  the  incubation 
of  the  artificial  substrate  samplers. 

Differences  in  reservoir  conditions  between  this  work  and  the  prior 
suiiiiier  1975  work  included:  1)  lower  water  temperatures,  and  2)  lower  reservoir 
stage.  This  latter  factor  siqni ficantly  restricted  mixina  in  the  upper  end 
of  Waconda  Bay  by  confining  the  flow  to  a relatively  narrow  channel  for  about 
3,000  feet. 

Midway  through  this  second  field  survey  the  Volunteer  Plant  was  closed 
by  the  Department  of  Defense.  Although  TNT  production  was  stopped,  wastes 
continued  to  be  released  into  Waconda  Bay  during  the  entire  period  of  study. 


FIGURE  1,  VICINITY  MAP  OF  VAAP  STUDY  AREA. 
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PHASE  I - FIELD  SURVEY 


The  purpose  of  the  December,  1976  field  survey  was:  1)  to  conduct  a 
dye  tracer  study  to  determine  flow  and  dispersion  patterns  in  Waconda  Bay; 

2)  to  make  on-site  measurements  of  D.O.,  temperature,  conductivity,  and  pH; 

3)  to  collect  samples  of  water  and  sediment  for  munitions  analysis;  and  4) 
to  obtain  limited  samples  of  benthic  and  periphyton  organisms  for  familiari- 
zation purposes.  Based  on  the  information  obtained,  a decision  could  be 
made  as  to  whether  or  not  a sufficient  concentration  of  munitions  existed  to 
justify  a more  intensive  biologic  survey,  and  if  so,  further  determine  the 
sampling  sites  locations  in  order  to  encounter  the  full  range  of  munitions 
concentrations.  The  survey  was  conducted  during  the  period  December  6-14, 

1976.  Daily  measurements  were  made  of  D.O.,  temperature,  conductivity,  and 
pH  for  the  5 day  period,  December  7-11  at  nine  sites  in  Waconda  Bay. 

Samples  for  munitions  analysis  were  taken  daily  from  these  nine  sites  plus 
the  effluent  from  the  last  treatment  pond.  Each  munitions  sample  consisted 
of  a composite  of  three  grabs  taken  during  the  day.  The  sites  utilized 

were  the  same  as  used  in  the  summer  1975  work  (see  Figure  2)  with  the  following 
exceptions.  Site  B-1  was  not  sampled  as  it  was  dry  land  at  the  low 

reservoir  levels  encountered.  One  sample  was  taken  in  the  middle  of  Transect 
F;  Transects  S,  T,  and  U were  not  sampled. 

Tabulations  of  the  field-measured  D.O.,  temperature,  and  conductivity 
values  are  shown  in  Appendix  Table  A-1.  Dissolved  oxygen  was  near  saturation 
values  at  all  stations  and  temperatures  were  in  the  5-lD°C  range.  High 
conductivity  values  were  observed  at  stations  near  the  waste  outfall  and  this 
parameter  appeared  to  be  a good  indicator  of  munitions  waste.  Equipment 
failure  precluded  measurement  of  pH.  Munitions  concentrations  found  in  water 
and  sediment  samples  are  tabulated  in  Tables  1 and  2.  The  concentrations 
found  exceeded  those  which  in  prior  field  studies  have  been  shown  to  produce 
effects  in  both  benthic  and  periphyton  communities  (Sullivan,  et  al . , 1977). 

The  dye  tracer  study  began  December  7 when  153  ml  of  Rhodamine  B was 
placed  in  the  plant  effluent  from  103D  to  1300  hours.  Based  on  previous 
experience  in  other  waters,  it  was  felt  that  a dye  concentration  of  near 
50  ppb  could  be  utilized  with  minimum  discoloration  of  the  lake  water.  Un- 
fortunately, the  lake  water  had  a light  green  tint  that  contrasted  sharply 
with  the  red  dye  making  concentrations  of  10-20  ppb  highly  visible.  State  of 
Tennessee  environmental  personnel  had  been  informed  and  approved  of  the  tracer 
study.  They  had  requested  however, that  obvious  discoloration  be  kept  to  a 
minimum  to  limit  citizen  complaint  and/or  alarm.  Dye  feed  began  at  about  32 
ppb  but  was  reduced  after  20  minutes  to  15  ppb  because  of  excess  color.  Even 
at  the  reduced  rate,  the  dye  plume  in  the  upper  end  of  Waconda  Bay  was  clearly 
visible  initially.  By  the  following  morning,  however,  maximum  dye  concentra- 
tion was  down  to  2.5  ppb,  and  was  not  visually  detectable.  By  mid-afternoon, 
maximum  concentration  was  2.0  ppb  and  most  of  the  dye  was  within  the  narrow 
channel  at  the  upper  end  of  the  bay.  The  rapid  decrease  in  dye  concentration 
observed  in  one  day  suggested  that  the  concentration  would  be  below  detection 
limits  when  it  moved  further  down  into  a much  wider  portion  of  the  bay.  Since 
it  was  desirable  to  observe  the  dye  further  down  bay,  an  additional  116  ml  of 
dye  was  released  into  the  center  of  the  plume  at  1600  hours  on  December  8. 

Even  with  this  addition,  the  maximum  concentration  found  the  following  morning 
was  2.0  ppb  and  the  peak  concentration  was  upbay  of  where  it  had  been  the 
previous  afternoon.  Transects  were  made  of  the  dye  plume  vertically  and 
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TABLE  1 

WATER  MUNITIONS  DATA.  MEAN  VALUES 
DECEMBER  1976 


HBHHIH 

VALUES  REPORTED  IN  uq/1 

TOTAL 

msfmsm 

TNT 

1 .3-DNB 

1.3.5-TNB 

EFFLUENT 

93 

53 

47 

4 

<1 

198 

A 

173 

89 

128 

6 

<1 

397 

B 

156 

88 

no 

5 

<1 

360 

C-1 

23 

18 

5 

<1 

<1 

48 

C-2 

78 

40 

34 

6 

<1 

159 

D-1 

9 

4 

<1 

1 

2 

17 

D-2 

14 

11 

<1 

<1 

2 

29 

E-1 

6 

3 

<1 

<1 

<1 

E-2 

2 

<1 

<1 

<1 

1 

_<  6 

F 

2 

<1 

<1 

<1 

<1 

< 6 

SEDIMENT  MUNITIONS  ANALYSES,  DECEMBER  1976 


horizontally  using  a continuously  measuring  and  recording  fluorometer.  The 
general  pattern  of  dye  movement  can  be  seen  in  Figure  A-1 . 

In  an  effort  to  better  understand  the  dye  movement,  reservoir  stage 
data  were  obtained  from  TVA  and  isoconcentration  plots  were  prepared  of  the 
conductivity  and  total  munitions  data.  These  data  are  shown  in  Figures  A-2 
toA-4.  Examination  of  stage  data  shows  that  the  dye  was  initially  introduced 
at  the  end  of  a rapid  drawdown  period.  From  December  7 through  December  10, 
the  reservoir  stage  increased  by  1.15  feet.  From  Transect  E upstream  this 
increase  represents  a volume  change  of  about  43  million  gallons.  This  is 
more  than  the  total  effluent  flow  from  VAAP  during  this  period  so  net  flow 
was  actually  upbay  at  Transect  E.  This  would  have  the  effect  of  holding 
effluent  in  the  upper  end  of  the  bay.  The  reservoir  level  dropped  approxi- 
mately 0.9  feet  during  the  morning  of  December  11,  which  had  the  effect  of 
moving  contaminants  downbay  at  an  accelerated  rate  and  which  the  conductivity 
and  total  munitions  data  demonstrate. 

Analysis  of  the  combined  data  shows  that  concentrations  of  pollutants 
in  the  first  3,000  feet  of  the  bay  were  roughly  0. 5-1.0  times  effluent  con- 
centrations. This  appears  to  be  due  to  the  limited  cross-section  of  the 
channel  at  the  low  reservoir  stages  seen  during  the  winter  months  (channel 
width  100  ft.  depth  5-6  ft.)  which  minimizes  mixing  with  cleaner  water  in  the 
bay.  Concentrations  then  drop  off  very  rapidly  over  the  next  1600  ft.  to 
levels  around  1-5  percent  of  effluent.  Variations  in  absolute  concentrations 
can  be  significantly  affected  by  changes  in  effluent  concentrations,  flow,  and 
by  changes  in  reservoir  stage. 

Four  sampling  sites  were  selected  for  "high"  concentration  ranges  in 
the  upper  2500  feet  of  the  bay;  three  sites  were  selected  in  "middle  to  low" 
ranges  within  1500  feet  of  the  point  where  the  bay  significantly  widens  and; 
two  additional  sites  were  selected  further  downbay  for  the  extreme  low  level 
range.  Three  sites  were  chosen  as  reference  points  in  the  adjacent  bay.  All 
sites  selected  for  the  second  phase  study  are  shown  in  Figure  3 . 
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Introduction 

During  the  period  March  1-5,  1977,  a second  field  survey  was  conducted 
in  Waconda  Bay.  Field  data  were  taken  daily  for  D.O.,  temperature,  and 
conductivity.  Dissolved  oxygen  was  near  saturation  values  at  all  stations 
and  temperatures  we  re  in  the  7-13‘’C  range.  High  conductivity  values  were 
observed  at  stations  near  the  waste  outfall  and  this  parameter  appeared  to 
be  a good  indicator  of  munitions  waste.  Equipment  failure  precluded  the 
measurement  of  pH.  Water  samples  were  taken  at  each  site  three  tinies  per 
day  and  composited.  Analyses  were  made  for  ammonia-N,  total  Kjeldahl  nitrogen, 
nitrate-N,  chloride,  total  hardness,  sulfate,  and  munitions.  Additionally, 
effluent  samples  (24  hour  composites)  were  taken  daily  during  March  for 
munitions  analysis. 

Sediment  samples  were  taken  at  each  site  in  Waconda  Bay  (three  samples 
per  site  where  each  sample  consisted  of  three  dredge  grabs).  These  samples 
were  analyzed  for  total  and  volatile  solids,  nitrate-N,  total  Kjeldahl 
nitrogen,  and  munitions. 

Except  for  the  munitions  samples,  all  analyses  were  performed  in 
accordance  with  procedures  of  Standards  Methods  (APHA,  1975),  (EPA,  1974),  or 
Chemistry  Laboratory  Manual  Bottom  Sediments  TTPA,  1969).  Details  of 
analytical  procedures  are  presented  in  Appendix  B. 

Characterization  of  Water  Qual i ty 


Mean  values  for  water  quality  data  are  shown  in  Table  3.  Typical 
concentration-distance  profiles  are  shown  in  Figures  4 through  6 for  nitrate- 
N,  ammonia-N,  and  chloride.  For  non-munitions  parameters,  these  data  show 
nearly  constant  average  concentrations  for  stations  1 through  5.  Sampling 
sites  1 through  4 are  located  in  the  relatively  narrow  and  shallow  channel 
existing  in  the  upper  end  of  Waconda  Bay  during  low  reservoir  stages. 

Station  5 is  at  the  iix)uth  of  this  channel.  The  channel  geometry  in  this 
area  limits  mixing  and  dispersion  of  effluent  with  other  water  in  the 
reservoir.  Down  bay  from  station  5 the  reservoir  widens  considerably  and 
concentrations  decrease  rapidly. 

Nitrate  nitrogen  levels  were  about  11  mg/1  at  stations  1 through  5 
dropping  to  1 mg/1  at  station  9.  Ammonia,  total  Kjeldahl  nitrogen,  chloride, 
total  hardness,  and  sulfate  also  showed  elevated  levels  at  stations  1 through 
5 with  concentrations  decreasing  from  stations  5 to  9. 

Five  specific  munitions  related  compounds  were  analyzed:  2,4,6  TNT; 

2,4  DNT,  2,6  DNT;  1,3  DNB,  and  1,3,5  TNB.  The  sunmation  of  these  five 
compounds  versus  distance  is  shown  in  Figure  7.  Mean  values  for  these  samples 
are  presented  in  Table  4.  Only  three  of  the  individual  lake  samples  showed 
the  presence  of  2,4  DNT.  Five  samples  contained  1,3  DNB.  This  compound  was 
not  detected  at  stations  6 through  9.  Most  samples  contained  2,6  DNT  and 
2,4,6  TNT  and  the  concentration  appeared  to  decrease  from  station  5 through 
9.  By  contrast,  1,3,5  TNB  was  more  often  found  at  down  bay  locations.  On 
only  two  occasions  was  this  compound  detected  at  stations  1 through  5 
compared  to  thirteen  occurrences  at  stations  6 through  9.  It  appears  that 
1,3,5  TNB  is  a breakdown  product,  probably  of  2,4,6  TNT.  Plant  effluent 
samples  contained  no  detectable  1,3,5  TNB  or  1.3  DNB. 


17 


STATION 

Mean  Values 

NUMBER 

NH3-N 

TKN 

NO3-N  Cl 

T.  Hard 

SO4 

mu 

1.23 

1 .91 

11.2 

17.2 

156 

155 

1.21 

2.01 

11.3 

16.9 

157 

162 

1.11 

1.84 

11.5 

16.3 

161 

157 

1.09 

1 .87 

11.3 

15.9 

159 

163 

5 

0.95 

1.59 

11.1 

15.2 

161 

206 

6 

0.44 

0.94 

5.59 

12.3 

125 

85.5 

mu 

0.23 

0.68 

3.24 

9.9 

107 

69.7 

8 

0.14 

0.50 

1.89 

8.7 

92.0 

44.2 

9 

0.07 

0.39 

1.04 

7.8 

82.0 

21.9 

EFF 

1.27 

1.97 

14.1 

18.3 

162 

164 

mgN/1 

mgN/1 

mgN/1 

mgCl/1 

mg/1  as 

1119504/ 1 

CaC03 

FIGURE  4.  NITRATE  NITROGEN,  MARCH  1-5,  1977 
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FIGURE  5.  AMMONIA  NITROGEN,  MARCH  1-5,  1977 
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TABLE  4 


MEAN  VALUES  FOR  MUNITIONS  ANALYSES  OF  DAILY 
SAMPLING  DURING  THE  PERIOD  MARCH  1-5,  1977 


STATION 

MEAN  Values 

in  u£/l 

NUMBER 

1 ,3-DNB 

1 ,3.5-TNB 

2,4-DNT 

2.6-DNT 

2.4,6-TNT 

TOTAL 

1 

4.1 

2.9 

< 0.10 

17.9 

48.0 

73.0 

2 

<0.25 

"0.75 

22.1 

38.7 

40.3 

102 

3 

<0.25 

vO.75 

V 0.10 

32.0 

43.7 

4 

0.4 

vO.75 

17.8 

37.2 

13.3 

69.5 

5 

0.5 

17.1 

V 0.10 

25.7 

26.2 

69.7 

6 

<0.25 

1 .2 

V 0.10 

12.1 

1.8 

J5.6 

7 

vO.25 

1 .6 

< 0.10 

7.6 

0.2 

::  0.8 

8 

<0.25 

10.6 

< 0.10 

1.3 

0.1 

:^12.4 

9 

<0.25 

20.6 

1 .3 

2.3 

0.3 

24.8 

24  hr  comp 
EFF 

<0.25 

/ 

o 

CJT 

28.3 

15.5 

44.5 

89.4 

stance 


FIGURE  7.  TOTAL  MUNITIONS,  f^RCH  1-5,  1977 


A station-to-station  comparison  was  made  using  the  following  chemical 
data:  NH3-N,  TKN,  NO3-N,  Cl,  total  hardness,  SO^,  and  total  munitions.  The 
data  for  each  parameter  was  first  converted  by  dividing  all  values  by  the 
maximum  value  in  the  set.  This  puts  between  a zero-to-one  range. 
Station-to-station  comparisons  were  made  by  calculating  an  average  Euclidean 
distance  as  follows: 


'AB 


s 


where 

n 


the  Euclidean  distance  between  Stations  A and  B 

the  number  of  chemical  or  chemical/biologic  parameters 
considered 

the  magnitude  of  the  i^  th  parameter  at  Station  A 
the  magnitude  of  the  i_  th  parameter  at  Station  B 


The  magnitude  of  the  Euclidean  distance,  A/^3,  increases  for  any  pair  of 
stations  as  the  number  of  parameters  considered  is  increased.  To  eliminate 
this  dependence,  an  average  distance,  d/\n,  may  be  calculated  (Sokal  and  Sneath, 
1963): 


AB 


where  d^g  = average  distance  between  Stations  A and  B 

A^g  = Euclidean  distance  between  Stations  A and  B 

n = number  of  chemical  or  chemi cal/biologi c parameters 
considered 


The  average  Euclidean  distance  values  were  used  to  construct  a phono- 
gram showing  the  relationships  between  the  nine  stations  (see  Appendix  for 
details).  This  phonogram.  Figures,  shows  stations  1 through  5 as  decidedly 
different  from  stations  6 through  9.  The  cophenetic  correlation  coefficient  of 
0.951  indicates  that  the  phonogram  is  an  accurate  representation  of  the 
relationship  between  stations.  Generally,  coefficients  greater  than  0.8  are 
considered  to  indicate  accurate  phonograms. 


Biologic  samplers  were  placed  in  Waconda  Bay  during  this  survey  and 
remained  for  a month.  During  this  period  samples  of  plant  effluent  were  taken 
by  VAAP  personnel  and  retained  for  subsequent  munitions  analysis.  This  was 
done  so  that  the  actual  munitions  concentrations  existing  during  the  incubation 
period  could  be  estimated.  The  results  are  shown  in  Figure  9 , and  Appendix 
Table  C-1. During  the  field  survey  period  (March  1-5)  the  effluent  munitions* 
averaged  89  pg/1.  During  the  following  six  days  effluent  munitions  concen- 
tration generally  increased  to  a maximum  of  about  2,400  ug/l  on  March  11.  It 
then  decreased  to  157  ug/1  on  March  17.  From  then  until  April  1,  average 
effluent  concentration  was  217  ug/l.  With  regard  to  the  specific  compounds, 

1,3  DNB  and  1,3,5  TNB  were  usually  below  detection  limits.  For  the  twenty 


* Total  of  the  five  specific  compounds  measured. 
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COEFFICIENT  0.951 


effluent  samples  taken  from  March  5 through  April  1,  1977  the  concentrations 
of  2,4  DNT,  2,6  DNT,  and  2,4,6  TNT  averaged  323,  98,  and  151  ug/1  respectively. 
Effluent  flow  rates  remained  relatively  constant  throughout  this  period  averaging 
5.9  MGD  (see  Appendix  Table  C-2). 

Consequently,  it  appears  that  average  concentrations  during  the  incu- 
bation period  exceeded  those  during  the  field  studies.  Total  munitions  in 
the  channel  extending  down  bay  to  and  including  station  5 probably  approached 
2,000  ug/1  and  averaged  near  500  ug/1.  Based  on  the  relative  munitions 
concentrations  found  at  station  6 through  9 during  the  early  March  survey, 
total  concentrations  of  munitions  may  have  risen  to  50  - 100  ug/1  at  these 
sites  later  in  March. 

Thus,  the  munitions  concentrations  encountered  by  the  organisms  colo- 
nizing the  artificial  substrates  may  have  been  significantly  higher  than  the 
concentrations  encountered  by  the  natural  substrate  organisms  prior  to  the 
early  March  sampling.  No  detailed  munitions  analyses  are  available  for  the 
period  proceeding  the  March  sampling. 

Characterization  of  Sediment  Quality 

Mean  values  for  sediment  quality  are  shown  in  Table  5.  Results  for 
nitrate,  total  kjeldahl  nitrogen,  and  total  solids  reveal  little  except 
slightly  lower  TKN  values  at  the  first  three  stations.  The  munitions  results 
show  that  1,3,5  TNB  and  2,4,6  TNT  clearly  predominate  1,3  DNB,  2,4  DNT,  and 
2,6  DNT.  The  former  compounds  were  found  at  all  stations  and  in  all  except 
one  replicate.  For  both  compounds  concentrations  at  stations  6 through  9 
were  higher  than  at  stations  1 through  5;  260  and  109  ug/kqfor  1,3,5  TNB 
and  97  and  133  uq/kgfor  2,4,6  TNT.  At  stations  1 and  2 only,  2,4  DNT  was 
found  at  <10  ug/kg  at  stations  1 through  5 only,  1,3  DNB  and  2,6  DNT  were 
found  at  --20  ug/kg. 


FIGURE  9.  TOTAL  rfJMTIONS  CONCENTPJi>TION  IN  VAAP  EFFLUENT  DURING  MARCH,  1977 


TABLE  5 


SEDIMENT  ANALYSES,  MARCH,  1977  j 

j 
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STATION 

NUMBER 

NO?- 

NO;, 

TKN 

TS 

1,3-DNB  1 

,3,5-TNB  2, 

,4-DNT 

2,6-DNT  2 

,4,6-TNT 

1 

<2 

350 

14 

<13 

73 

n 

10 

91 

2 

4 

333 

12 

<12 

<179 

<7.9 

17 

119 

£2 

352 

17 

<14 

76 
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PHASE  11  PERIPHYTON 


Introduction 

The  periphyton  or  "Aufwuchs"  coninunity  is  an  assemblage  of  attached 
microorganisms  (primarily  algae)  growing  on  the  surface  of  submerged 
substrates  and  forming  a coating  which  is  coimionly  a green  or  brown  color. 

This  coninunity  consists  of  both  autotrophic  (i.e.  unicellular  or  filamentous 
algae)  and  heterotrophic  (bacteria,  protozoa,  rotifers,  ect.)  forms. 

These  organisms  are  well  suited  to  biological  investigations  since 
they  remain  at  fixed  lo.-ations  and  are  sensitive  to  environmental  altera- 
tions, Their  populations  and  biomass  are  relatively  easy  to  quantify  using 
standard  laboratory  procedures  and  are  adaptable  to  a variety  of  statistical 
analyses.  The  results  of  the  survey  conducted  at  VAAP  during  the  winter  of 
1977  consist  of  two  parts:  1)  data  pertaining  to  the  field  survey  and  2) 
treatment  of  this  information  via  statistical  inference. 

The  community  structure  was  examined  using  diatoms  and  filamentous 
algae,  chlorophyll  biomass  and  vital  staining.  Use  of  the  staining  tech- 
nique showed  the  live-dead  relationship  among  the  periphyton  population. 

The  number  of  diatom  frustules  that  are  inactive  in  a coninunity  colonizing 
a glass  slide  cannot  be  estimated  using  standard  techniques.  Since  the 
inactive  portion  may  be  substantial, an  estimate  was  made  on  selected  stations 
using  tetrazolium  violet.  The  methodology  for  the  various  procedures  is 
appended. 

Artificial  Substrate  Colonization  Studies 

Glass  microscope  slides  with  a surface  area  of  3871  iiin^  were  exposed 
and  allowed  to  colonize  for  30  days  at  the  12  selected  stations.  Periphyto- 
meters  containing  8 slide-unit-cartridges  were  placed  together  in  suspension 
units  as  shown  in  Appendix  D.  Suspension  units  containing  both  6 and  8 
peri  phytometer  units  were  used  in  the  study.  The  suspension  units  with  the 
lower  number  of  peri  phytometers  were  used  to  provide  additional  slides  for 
the  statistical  portion  of  the  investigation  and  were  added  at  those  stations 
predicted  to  reflect  high,  medium  and  low  stress  from  VAAP  effluent. 

Organism  density  ranged  from  27  per  nin^  at  Station  2 to  37,267  at  Station 
11.  Numbers  of  species  per  station  ranged  from  22  to  52  in  Waconda  Bay  and 
from  76  to  83  in  the  reference  bay.  To  simplify  data  presentation  and  elimin- 
ate rare  occurring  species,  a culling  routine  was  completed  and  all  species 
whose  numbers  did  not  constitute  at  least  one  percent  of  the  total  population 
at  any  one  station  were  eliminated.  This  procedure  resulted  in  a reduction  of 
species  from  124  to  30.  However,  the  maximum  drop  in  population  density  at 
any  one  station  was  7 percent.  Table  6 shows  the  differences  in  density 
and  number  of  species  before  and  after  culling.  The  culled  data  are  shown 
in  Tabi<  7.  Organism  density  is  significantly  different  at  stations  1 
through  5 as  compared  to  stations  6 through  12  (see  statistical  report  in 
Appendix) . 


29 


TABLE  6 


A COMPARISON  OF  POPULATION  SIZE  AND  NUMBERS  OF  SPECIES 
IN  CULLED  AND  UNCULLED  DATA  AT  U. 

PERIPHYTON  - ARTIFICIAL  SUBSTRATE 


STATION  1 


CULLED 


no/niir 


UNCULLED 
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no/niiT 
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Z DIFFERENCE 
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no/nm 
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The  dominant  species  in  Waconda  Bay  are: 

Achnanthes  minutissima 
CjniibeIJa  prostrata 
Diatoma  tenue  var.  e Tonga turn 
Fragi laria  vaucheriae 
Helosira  varians 
Synedra  rumpens 

These  species  reached  population  levels  of  at  least  1,000  individuals/ 
mm^  at  one  or  more  stations  in  Waconda  Bay.  Diatom  populations  were  reduced  in 
the  impact  area  from  stations  1 to  5 although  the  actual  number  of  species 
was  approximately  the  same  as  down  bay. 

The  marked  population  increase  at  station  6 can  be  correlated  with  the 
water  chemistry,  which,  during  the  March  1-5  survey,  showed  a significant  drop 
in  concentration  of  several  components  between  stations  5 and  6.  (See 
chemistry  section.)  Considering  the  dominants,  the  ratio  of  density  for 
dominant  organisms  between  stations  6 and  5 where  the  Bay  widens  is  as  follows: 


Achnanthes  minutissima 

- 48 

Diatoma  tenue 

- 25 

Fraqilaria  vaucheriae 

- 14 

Melosira  varians 

- 9 

Synedra  rumpens 

- 51 

Figure  10  shows  the  clustering  of  stations  by  the  Pinkham-Pearson  Index 
using  culled  data  with  mutual  absence  important.  The  phenogram  shows  that 
stations  1 through  5 are  decidedly  dissimilar  from  stations  6 through  12. 
Additionally,  the  reference  bay  stations,  11  and  12,  form  a sub-group  somewhat 
different  from  Waconda  Bay  stations  6 through  9. 

Shannon-Weaver  diversity  values  are  shown  in  Table  8.  Theoretically,  the 
higher  the  diversity  the  "healthier"  the  community.  However,  examination  of 
the  results  shows  that  diversity  in  the  area  of  maximum  stress,  stations  1 to 
5,  is  higher  than  at  the  downbay  stations  6 to  9.  This  anomolous  condition  is 
caused  by  the  severe  drop  in  population  density  of  all  dominant  species  at 
stations  1 to  5 which  results  in  a condition  of  better  balance  between  the 
surviving  species.  Since  diversity  considers  only  relative  proportions  of 
species  and  ignores  absolute  population  density,  it  seems  to  be  of  limited  value 
in  a situation  where  the  main  input  is  seen  as  reduction  in  population  density. 

Natural  SubstraJ.e  Col onijation  StjcH^s 

Natural  substrate  collections  were  made  at  9 stations  in  Waconda  Bay 
and  3 stations  in  the  reference  bay.  A description  of  these  samples  and  the 


32 


PAGE  33  - FIGURE  10  Is  also  PAGE  167  - FIGURE  G-1 
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Pinkham-Pearson  Biotic  Similarity 


COPHEflETIC  CORRELATION  COEFFICIENT  .871. 


SHANNON-WEAVER  DIVERSITY  INDICES  (H)  FOR  VAAP  ARTIFICIAL  SUBSTRATE  DIATOMS 

MARCH,  1977 


resulting  data  are  shown  in  Table  9 and  10.  Generally  scrapings  were  made 
from  submerged  branches,  rocks,  and  styrofoam  floats.  Three  substrates  were 
sampled  at  each  location  except  stations  4 and  8.  Counts  were  made  to  approx- 
imate 500  diatom  valves.  Comparisons  among  stations  can  only  be  made  relative 
to  numbers  of  individual  species  colonizing  these  substrates.  The  overall 
total  number  of  species  ranged  from  11  to  77  with  no  marked  differences 
related  to  station  locations.  In  fact,  the  two  extreme  values  both  occurred 
at  station  3.  The  data  are  culled  to  eliminate  species  that  don't  make  up  at 
least  one  percent  of  the  total  at  at  least  one  station. 

The  five  most  commonly  occurring  species  were: 

1 ) Achnanthes  minutissima 

2)  Fragilaria  capucina 

3)  Fragilaria  vaucheriae 

4)  Stephanodiscus  invisitatus  ■ 

5)  Synedra  rumpens 

In  considering  the  distribution  of  individual  species  it  is  interesting  to 
note  that  Achnanthes  minutissima  was  as  common  an  occurence  in  the  impact 
zone  as  in  other  areas  of  Waconda  Bay  and  the  reference  bay.  This  suggests 
that  Achnanthes  is  not  as  susceptible  to  munitions  waste  as  might  be  expected. 

In  an  earlier  study  conducted  during  the  summer  of  1975, Achnanthes  minutissima 
showed  a similar  response  and  accounted  for  80  percent  of  the  diatom  population 
colonizing  artificial  substrates  in  areas  of  severe  stress  from  VAAP  waste. 

The  two  Fragi laria  sp.  reached  a maximum  of  approximately  one  fifth  the 
population.  These  levels  were  attained  only  at  downbay  stations.  With  the 
exception  of  Synedra  rumpens  at  station  4, the  minor  dominants  were  severely 
growth  restricted  in  the  high  stress  zone  between  stations  1 through  5. 

Pinkham-Pearson  comparisons  among  the  various  samples  are  shown  in  Figure 
11.  Clustering  of  stations  using  natural  substrates  results  in  a grouping 
based  primarily  on  substrate  type  rather  than  location  relative  to  the  waste 
discharge.  For  example,  the  first  station  cluster,  with  one  exception,  has  a 
common  denominator  of  wood  scrapings  as  a substrate  type.  The  second,  begin- 
ning with  4AR  includes  predominately  rock  scrapings.  Styrofoam  substrates 
principally  compose  the  next  group  and  include  samples  from  stations  6 through 
8,  11  and  12,  Treatment  of  the  data  by  the  Pinkham-Pearson  analysis  is 
especially  useful  in  this  case  to  point  out  the  substrate  dependence  of  peri- 
phyton diatoms.  The  phenographic  display  in  this  instance  further  demonstrates 
the  utility  of  uniform  substrates  for  periphyton  investigations.  The  cophenetic 
correlation  coefficient  of  0.734,  however,  indicates  that  caution  should 
be  exercised  in  interpreting  this  phenogram. 

Filamentous  Algae 

Filamentous  forms  colonizing  glass  slides  were  restricted  to  green  and 
blue-green  species.  The  data  are  shown  in  Table  11.  Six  genera  predominated 
with  the  blue-green,  Schizothrix  calcicola,  the  most  abundant  organism  at  all 
stations.  The  second  dominant,  Osci 1 1 atoria  sp.,  was  observed  at  all  sampling 
locations  but  occurred  only  sporadically  between  stations  1 through  5.  An 
Anabaena  species  was  recorded  only  once  from  station  12  in  the  reference  area. 
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TABLE  9 


INVENTORY  OF  NATURAL  PERIPHYTON  SUBSTRATE 
SAMPLES  COLLECTED  FEBRUARY  28,  1977 


Station 

Substrate 

Description 

Waste  Pond  Spillway 

Rock  Scraping 

Waste  Outfall  Stream 

Rock  Scraping 

"Old  Station  A"* 

Branch  Scraping 

"Old  Station  A" 

Branch  Scraping 

"Old  Station  A" 

Branch  Scraping 

S-1 

Branch  Scraping'^ 

S-1 

Branch  Scraping^ 

S-1 

Log  ScrapingW 

S-2 

Branch  Scraping^ 

S-2 

Branch  Scraping'^ 

S-2 

Branch  Scraping'^ 

S-3 

Rock  Scraping*^ 

S-3 

Branch  Scraping*"^ 

S-3 

p 

Rock  Scraping 

S-4** 

Dock  Scraping^'* 

S-5 

p 

Rock  Scraping 

S-5 

Rock  Scraping^ 

S-5 

P 

Rock  Scraping 

S-6 

U 

Stump  Scraping 

S-6 

Rock  Scraping*^ 

S-6 

Dock  Styrofoam^ 

Substrate 

Station  Description 


S-7 

c 

Dock  Styrofoam'^ 

S-7 

Dock  Styrofoam^ 

S-7 

Branch  Scraping'^ 

S-7 

Stump  Scraping'^ 

S-8 

Styrofoam  from 
"No  Wake"  Marker^ 

S-8 

Scraping  from 
Metal  Buoy^' 

S-9 

Scraping  of 

Metal  Buoy^' 

S-9 

Rock  Scraping*^ 

S-9 

p 

Rock  Scraping'' 

S-10 

c 

Dock  Styrofoam 

S-10 

Rock  + Old  Bottle 
Scraping*^ 

S-10 

Rock  Scraping*^ 

S-10 

Surface  Sediment*^ 

S-11 

Wood  Dock^ 

S-11 

Dock  Styrofoam^ 

S-11 

Rock  Scraping^ 

S-12 

Rock  Scraping*^ 

S-12 

C 

Dock  Styrofoanr 

S-12 

Branch  Scraping'** 

*20  yards  from  VAAP  outfall 

**Only  1 substrate  available  at  Station  S-4 


W = Wood;  R = Rock;  S = Styrofoam;  M = Metal;  B = Surface  Sediment 
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Taxonomic  Classification 


B»CUL»HIOI>MlfT»  (OUTOhS) 


iCW»-»NT«E3  t‘I6U*  V HtTEIQVtLve 
iC^NJ'iTHES  L*'»CtOL»T»  V OUBIi 
HINUTISSIM* 

»CM:,Is7Mt3  SOLLII 
jKor'iif  ONf  IS 
SSTfilOnHl.*  FUKmQS* 

e»L0NtI.S  MTiLlN* 

CDCCOWEIS  ►'i.iCtuTbLA  v EUGLr'TA 
CVCLOTeuLA  »IQ"U3 

CYCLCnELL*  MtNEGKiMiNi 
Ct^'detLA  OtLiCATuLA 

cv^eebL*  bic^ictPhAL* 

CV^BfLt*  Pr»0«T»9ATA  V AUEPSh»AlDII 

CVHmFULA  tomida 
OIaTJma  TtNut  V ELONGaTum 

fu^aTIA  PtCTiNALIS  V hTNOH 
FRaGIlARIA  CAPuCINA  V hESOLtPTA 


cvhmfula  Tomida 
oiaTjma  TtNut  V eloncaTum 

fu^aTIA  PECTINaLIS  V mtnoh 

wtiili  '' 

FPAGTLAWIA  C'J^SIRUF.NS  Ki 
FfiA^J^AMjA  ^UJ*Sr«UtNS  V ventew 

FRAGTLAklA  VAUCPeWlAt 
fiOMPMQNfMA  aNGuSTATuK 
^O^'PMO'^t^A  PAMVOLUM 

mFLQSiRa  OISTanS 
•*funsJ»A  J^A'^VW*^* 
mELUSIRa  Iyalica 

NfcLaSIKA  VAWlAkS 

maVICULa  CAPIIATA  V HONQARICA 

MAVICULA  CInCTA 

NaVICULa  CUnFERACEA  V PtPEGHiNA 
NAVJCULA  COnTLnTA  V »1CEPS 
MAVfci^lA  C'<y^lOCfcPH*LA 

MAVIpULa  HtNiSCOLUS  V UPSACIENMS 
NaVICOLA  minima 
NAVICULa  HUTICA 

MAVKULA  PELLICULQS* 

NAVirilLA  PoPUUA  V MOTATA 
NAVICUCA  NAOIUSA  N 

Klvi^uti  ShvncScey^mIEa^v^GEHhani 

NAVICOLA  5A1.1NAHXUM  V INTERMEDIA 

NAVICULa  StCMfcTA  V aPICuLATA 
NAVJCJLA  3tMlNULUM 
NAVICULa  SIMULA 

NAVICULa  5U8HAMUI  aTA 

mmiut 

NAVIC^l-*  SP  A 
NTTMcmIA  ArCUMOOATA 
NIT23CN1A  AcCUMUOATA 


N T7APMJA  AHMHIoIA 
N tJsCxi*  HTr» 

N Tjscwi*  OISSIP»T» 


NITISCMI*  Ot'l’ICUL* 

NlT/SCXl*  f ItI»'3S»1IS 
NITlscwl*  FRUSTULU**  V 8U»S»1.1R» 

RITI9CMJ*  KljtTI 
RIIJSCRl*  RlLt* 

NiTllCHI*  RtCT* 


nJtJuCmIa  aA  b 
sTCPMANOOIbc'JS  ASTNaea 

STCPmaNOQISCUS  aSTNAEA  y,  MlN 
STkPMANuSlscUj  iNVISlTAtUS 
SunXnELLA  U«aTa 

SUNiPtLLA  SutCiCA 
SVNfoMA  OELtCAtlSSlMA 
•VNtOHA  BUM^tNS 

SyNfORA  RuMpENS  y FAMXlXAhIE 
SYNenRA  Ulna 
TAeiLlARU  PLUCCULObA 


total  numbeH  up  CRGanIS^a 

NUMpEH  UP  TaVA 


TABLE  10.  (Continued) 


Number  of  Organisms  at  station: 


Taxonomic  Classification 


Number  of  Organisms  at  station: 

• ?»•  • fll  • Ten  • • »»»•  * 


*C^NastHt3  V **t  )'•  *l 

THt»  ^ OiJHl* 

♦ C^^NAsTwtj  MiSuilbM*’4 

AS'.'H  H Ofxt  1 a * 1 1 Nt  * 

CALONtlo 

Cl^rCONt  IS  »»i.4Ctsr,,L4  » fcl>*»LT»»T4 
cvri.on.tL4  4^o^os 

CVCtOTttt.4  Ht  \l 
f yqton cl*  >ULLr»>''4 
C ti  * *► » i • I ? 

Cv»'j.t  1 1 * 'Jr  t U 4 f\.rc4 
CV*’1*1  k.  L A w *}  V I '' 

Cv»«!»rtLA 


rVMr*ttt*  ►‘xoS'WaTa  V 4iit  *»S»  *tO  I J 

CV*<M^tL4  »UMlU4 

'MAli‘»4  TtNjt  V tLON-.aTy 


rU'^iTl*  ►'tCT  l '<4t  IS  » «InON 
4«*yfL*-<IA  C*'0C!'‘4  V «t<JLl**T4 

If  It  *N  I * c*^uc  i •* 
rKA..  tc  A'' I 4 CONMNir'*'  •• 

M*v.it*f'L*  ruN8ti4ttN:4  V VtSTfN 
I*  H lb  tt  *>«  I 4 *1  ^N4  r A 

►••J  41.  K 4 -<  ! 4 V *'JCH»  L 1 At 
A'*40iTAtu’4 
OijMf'u.lNt”*  P*'4VUtOM 


’<tLO»IN*  y»Hl»Ni, 

NiXICOl.*  ClPl'Af*  V «1/’>6*«ICA 
Xi«!CUl.t  CINCU 

NiVIC'lL*  CtiN'tfACl*  V •'f“£r.HlN» 
CiInUM*  V mc|l’4 

NAVlfULA  •'tNlaouLUS  * UlASALlfNtlS 
V*»iCUl-A  MINI*'* 

N»VlCuLi  "UTU* 

NAVICULA  PtLl.lCul.USA 
navIcuLa  pupula  V puTaia 
NAVKuLA  PAplJSA  N 


I CutA  |JA0IUSA  X ftlSttCA 
( tuLA  hpyNCUCtPMALA  V UtHPANl 
< CULA  SA^INaKIUM  V IMTtMntUlA 


NA*ICUI.A  PtC^tTA  V aPICULATA 


VAVirULA  SUbpaHULAIa 

NA^JCULA  }U,(' Ii,i5.tthA 

AAPffUlA  r"ftAlAABMI 


NAVICULA  SP  A, 

"'t’i'CPl*  ACUSpODATa 
NjTJStPlA  acCupuOaTa 

NIT/ACPIA  APPPlBIA 

snis?::!*.  ?liiip.TA 

UjTiACPIA  OAPtlCULA 

2l;i'56ni*.  fiti'TSi!5i\  SUM.LIN. 

Nlt/SCPIA  AjttllNOlAPA 
xlfisCPIA  PAllA 
n!*}sO"Ia  PtCiA 

PlflSrpIA  S^nuaTa  V TaPILUAHIA 

sTr^ptPuC<I.scp^  asTpaca 

ST*  PPAPUCp-"?  ASt»A*A  y.  PINUTUA. 
At*  PpaPuM  n I'A  INVISITATIIA 
SUPIPUUA  O^AlA 

AUPIPILLA  SiitCtCA 

V # 4»»lt.  I4KIS 

|TNfON4  IILN* 

rASiLLA^lA  fLUCCULUSA 


stf PpApuPI.'cpS  Ajt 


TuTAA  PUPAfP  UP  OPISAhlSPS 
PUMSIP  UP  TaAA 


TABLE  10- (Continued ) 
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Taxonomic  Classification 


tt*CIl,L4HiO#*HrTA  (OUTOHa) 


4C«Ni\T^»ea  LA^CfeOV 
»C»*N*sTHt3  fllNuTu 

ifNMl^ThfS  SOLLll 

4**p«ntO''«Cls  ''ItnC* 
43Tt**lPNf.Ll.*  fUhHO 


4 V HtTCROVALVC 
i3uAr4  V DUBXA 
USiHA 


CAUONf  Ji  rttALlNA 

fOCC^^tlS  *'L*^tNTUL4  V tUGLYPTA 
CYCLUttLL*  4TJ‘^U3 

CYCLOTfLlA  MtNtuHiKlANA 
C»CLi'rUL4  attLLlGfH4 
CrWbU.lA  i^lS 

CYHyttLA  UtilCATULA 
CV’^btLlA  LAtvIa 
CY'^bCLlA  --  - - 


a 

CLRHAUA 


CV^HFLLA  PHoaTRiTA  V AUgMSNALOlI 

CY'^BtLLA  TUMlOA 

OIATOma  TtNut  V ELQiyGaTuH 

KUNOtlA  PLCTI'NALIS  V HInOH 
^RAuILAPIA  CA^!^.'C^^A  V heaottPTA 
PHAGtLASlA  CA^'UCX^A 

PRAGTLANtA  N 

fR4&|L4-<I4  V Ve-<"€R 

rPAuTLAPlA  PINNaTA 

phaoIlakIa  vauchfRIaC 
GLlHPM?\t^A  A^bUaTATUR 
GO^>'hO%l'^A  PAlVu^_ijM 

^(LOSj-^a  t-lSfANS 
•*LLLiaI«A  G’Ya'OJLATA 
^ELOaiKA  XTaLICA 

MLLOSIP4  V4R1AN8 

NaVICuLA  CAPiTArA  V HUN6ARIC4 

RAVtCUlA  CInCM 

NaVICULa  COnFLRaCEA  V PEH{QRIN4 

utmt 

NAViroLA  H6N13CULU3  V UR3ALlPNaia 
RAVrCUCA  rtiNirtA 
NAVlCJlA  '’'JTICA 

NAYXCUlA  PFluICULOAA 
NAVXruLA  PuPULA  V MuTATA 
RAVlfi-LA  RAO  I US  A h 

NAVJf.UA  ■'^OI'JSA  V TtNfLLA 
NAVJC  LA  •*•«>  -tOtl  PMAt  A V utrtMANi 
NAVJCi.uA  ^AuINA^lgM  V iNTtMHfcOIA 

NAVir.:^.A  br^-’LTA  V APiCUlATA 
NAvirjLA  MN'jlu'‘ 

NAVfCULA  CF  AImClA 

NAVJCuLA  aUH"A«gt.ATA 
NAVjCi'LA  juuTjj  J3ajrtA 
NAVICULA  TPJtNCPANNil 

NAViriiiA  3P  h 
Njt^XCHlA  4rCJ‘'QpATA 
NIT/SC'^IA  ACCUNOSAtA 

N I T / AfN  I A 4H»'M  JiJ  I 4 
NlWU»*i4  rtjtA 
NlTXSgMlA  l'I^jIPaTa 


1 T;SCh  a OANTlQytA 
^ ^aSl-^  a ailILunnij 
w fJSfN  4 #fi'>STOLu*» 


V ajbaAitNA 


Nl  T ; sc  1 A A .It  T X j ng  I ana 
N?r*'SC«jA  -AiAA 
n|t^sc*<ia  -m  Ta 


NlT/STHtA 

Nff.sfilA  3P  ^ 
sTf  Pnan^t  : it  ''■» 


ATa  V tapillanta 


IffP-ANiUpr  ^ ASTnaFa  V, 
4TtP*«ANjM»C‘'f*  *NVta!TATu3 

lUPl^l  LI  A a>  A U 

ai’PlVlLL  A P.  ' V lu  A 
SYStPHA  ntul  aM8.U"A 
3YNt  ^11.  A N ,,  .rS  *5 

3YNtO*<l  M -ptNS  V AAMiLtA^lS 
SYnI^NA  l.g.<4 
TAttlLLAxlA  ALUCCULO^A 


fw^AL  N(iM0(Mi  gr  ONt.ANJSnS 
NUNfttN  UF  TaVA 


Number  of  Organisms  at  station: 
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Taxonomic  Classification 


C0l*^0H5) 


AC»<N*NThIS  l«UUA  V HtTf.NjvALVf 
V Oubl* 

4C‘^S*  T„t3  HlNOTl^SlHA 
Th»  S SOtL I I 

isnMfu  *u.t  13  vITh*.  a 
A5Tt»<I  JSk  lLA  f UMH05A 

TALON*  13  H»4t JV* 

CUCCTNUS  ►‘uALtNToLA  • tU&LVPTA 
CtCLOlfcuLA  4*uHij5 

CVCLOTH.LA  MbNLliHlNiA'^A 
CVCLf^tELLA  SULLlbfhA 
CVMBtLLA  AP^ INIS 

Cr**atLL4  OtLiLATuLA 
cr«5^LU*  L*evZ3 
CV^^BClL*  Mlc^UCtPHALA 

CV^tffLL*  pnoSThata  V AUCNSWALOII 

CY^JfLL#.  TumIDa 

OlATHhA  TtSut  V iUONGATuH 

eUNOTlA  PkCTlNALlb  V 
FRaGIlAnIa  CA^’UCl-A  V *^L&ULtPTA 
FflAGfLAHlA  C*POClNA 

FPAGTLANIa  C^^^T**UFnS  n 
► »A5nA^TA  rjNaT*^utN5  V VEntCH 
FRAurLANU  PIN^ATa 

FPaGTLAL'Ta  VAUCMtHlAE 
ROMPNONt^’A  ANu'JSTAfoM 

GOmPnUnl*’*  PANVUuub 

«|LQ2i«A  OI^tANS 
H|l5sNa  f.RANULATA 
*<lL59i><A  lULlCA 

»<ELOSJPa  vawUnj> 

NAViruLA  CAPMaTa  V MUNCAKICA 
NAViCULA  CInCYA 

NA»^ICULA  Ct>N»'iPACtA  V PFHtCPiNA 
navICJLa  CONtbNTA  V BICEPS 
NAVJCULA  C^YP'OCEPhAL* 

naVICULA  MtNldCoLUS  V UPSalIEKSIS 

NAVir^LA 

NAVICUIA  'iUTlCA 

NAVICUIA  PtiLlCOLOSA 
NAVTCVLA  PuPULA  V hutata 
NAVfCJLA  HAOIJSA  N 

navjcula  naoioSa  y tenella 

NAVICjLA  NHvNCUr  tPHALA  V (it'iMANl 
navICULa  SalInahiu'^  V INTERmeoia 


navicula  ^tc**uTA  y apiCulaya 

NAVICULA  StMtNLL'jP 
NAVfCuLA  CF  AZNJlA 


NAVICULA  SlirfNAMyLATA 
NAVJC  )LA  SUtf T u jSsiHA 
NAVjCULA  THltNEHANNlI 

NAVICoLA  SP  b . 


NAVICOLA  SP  b . 
"ITMCXI*  iMt-rtlBl 

sHl?s:iJ  SiiJip. 

NITIJCHI*  OfNTJrul 

:i\Wi 


Utl)'  V SoOSlLtNl 


Njrzscxj*  KutWixtu'n 


xl»  KjLti 
Xjj 


S^'IU»T*  V T«BELL>01t 

SuSixlCl*  uy»iA 

SU»I«I[LL*  SutLJC* 


STNECNt  Myxptxs  y EiXUH^IS 
SrxiOM*  III.N* 
riBCLLtxK  ykJCCULOK 


TuTtk  NOntitN  ur  ORStNlSNS 
x<J*1dCH  uR  Ttit 
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Number  of  Organist  s et  station: 


UiM  ♦ UCM  ♦ UAi| 


UBI  « XiCb  • 


I 


TABLE  11 

PRESENCE  - ABSENCE  DATA  FOR  FILAMENTOUS  ORGANISMS 
COLLECTED  FROM  VAAP  ARTIFICIAL  SUBSTRATES  - MARCH,  1977 


SPECIES 


Schizothrix  calcicola* 


ChloroDhyceae 


Mougeotia  sp  1 
Rhizoclonium  sp  1 
Stigeoclonium  spl 


* dominant  at  all  stations 


Three  species  of  Chlorophyta  were  present.  Of  these, only  Stjgeocloniuin 
occurred  in  the  impact  zone  between  stations  1 through  5.  Green  algae  repre- 
sented  by  Stigeoclonium,  Rhizocloniuni  and  Mougeotia  were  found  more  frequently 
in  the  reference  bay.  Only  Stigeoclonium  was  present  in  Waconda  Bay  samples. 

Chlorophyll  8 and  Organic  Biomass 

Estimates  of  plant  pigment  and  biomass  levels  were  completed  at  all 
stations  except  S-10,  where  the  peri  phytometer  suspension  unit  was  lost. 

Sixteen  replicate  slides  for  each  of  the  two  analyses  were  collected  from  the 
suspension  units.  These  data  are  presented  for  each  station  in  Figures  12 
and  13.  Table  12  contains  suiuiiary  information  for  means  and  ranges  as  well 
as  the  calculated  Autotrophic  Index  based  on  mean  values  per  station. 

Chlorophyll  a concentration  on  colonized  slides  at  the  various  stations  re- 
vealed a pattern  consistent  with  population  counts.  Stations  1 through  5 had 
virtually  no  chlorophyll  but  a sharp  increase  occurred  at  S-6.  A maximum 
average  concentration  of  AOmg/m^  was  observed  at  Station  7.  Pigment  levels 
gradually  declined  at  downbay  stations  to  17  mg/m^  at  S-9.  The  reference  bay 
samples  showed  means  of  14  and  28  and  could  be  considered  comparable  to  outer 
Waconda  Bay. 

Organic  biomass  shows  a similar  pattern.  Some  heterotrophic  growth 
occurs  in  the  high  stress  zone  between  stations  1 through  5.  Further  downbay 
periphyton  biomass  increases  with  a maximum  of  4.8  gms/nr  at  S-8.  A maximum 
mean  level  of  6.2  mg/m2  was  observed  at  S-11  in  the  Reference  Bay. 

Calculation  of  the  autotrophic  index  suggests  an  imbalanced  environment 
among  all  stations  except  S-7.  According  to  Weber  (1973),  ratios  greater  than 
100  indicate  organic  pollution.  Autotrophic  index  values  for  the  winter  survey 
ranged  from  90  to  55,000.  Maximum  ratios  were  observed  at  stations  1 through  3. 

Vital  Stain 

Three  replicate  slides  from  the  peri  phytometer  units  were  exposed  to 
tetrazolium  violet  to  estimate  the  physiological  activity  of  diatoms.  Cells 
were  separated  into  live,  senescent, or  dead, based  on  the  uptake  of  dye  and  the 
presence  of  chlorophyll.  Cells  completely  filled  with  stain  were  classified 
as  live.  Those  which  contained  chlorophyll  but  no  stain  were  considei'ed  senescent, 
Those  with  neither  dye  nor  chlorophyll  were  classified  as  dead.  The  technique 
offers  some  insight  into  growth  characteristics  of  diatoms.  Present  methods 
of  diatom  enumeration  using  cleaned  frustules  probably  leads  to  overestimates 
in  that  all  cells  are  considered  to  be  viable. 

The  data  have  been  expressed  as  fractions  using  each  category  as  it  is 
related  to  the  total  number  of  individuals  observed.  This  infotmtation  as  shown 
in  Table  13  reveals  that  conditions  existing  in  the  channel  between  stations 
1 and  4 resulted  in  nonviable  populations  that  ranged  up  to  46  percent.  A 
sharp  demarcation  occurred  at  station  5 where  the  percent  dead  organisms 
dropped  to  7 percent.  A further  reduction  to  1 percent  was  noted  at  station 
6 where  93  percent  of  the  organisms  were  vital  staining. 


44 


Down  Bay  Distance  (miles) 


FIGURE  12.  PERIPHYTON  CHLOROPHYLL  a,  ARTIFICIAL  SUBSTRATE,  MARCH  1977 


Down  Bay  Distance  (miles) 
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FIGURE  13.  PERIPHi'TON  BIOMASS,  ARTIFICIAL  SUBSTRATE  , MARCH  1977 


TABLE  12 

PERIPHYTON  CHLOROPHYLL  a,  BIOMASS  AND 
AUTOTROPHIC  INDEX,  MARCH  1977 


Station 

2 

Chlorophyll  a mo/m 

2 

Organic  Biomass  gm/m 

Autotrophic  Index 

Mean 

Range 

Mean 

Range 

1 

.02 

.01-. 06 

1.1 

.61-2.2 

55,000 

2 

.02 

.01-. 04 

.53 

.27-1.0 

26,500 

3 

.08 

.02-. 1 

.34 

.02-1.3 

4,250 

4 

■ .2 

.1  -.4 

.10 

.02-  .26 

500 

5 

.5 

.2  -.8 

.50 

.22-1.0 

1 ,000 

6 

14 

7 - 19 

2.3 

1.2  -2.8 

164 

7 

41 

29  - 52 

3.7 

bh 

90 

8 

23 

13  - 34 

4.8 

209 

9 

17 

5 - 35 

3.4 

.57-6.1 

200 

10 

- 

- 

- 

- 

11 

28 

16  - 44 

6.2 

4.7  -8.2 

221 

12 


14 


11  - 18 


-6.8 


336 


TABLE  13 


PHASE  II  - MACROINVERTEBRATES 


* 


Introduction 

Aquatic  macroinvertebrates  are  a diverse  group  of  small  aquatic  animals 
comprised  of  snails,  clams,  arthropods,  annelids  (segmented  worms  and  leeches), 
planarians,  and  coelenterates.  Of  these,  oligochaetes  and  chironomid  (midge 
fly)  larvae  account  for  the  majority  of  the  organisms  in  this  study. 

Aquatic  macroinvertebrates  are  a major  biological  component  of  aquatic 
systems  and  form  an  important  part  of  the  food  chain.  They  feed  on  detritus 
and  microscopic  plants  and  animals.  They,  in  turn,  are  eaten  by  small  fish 
which  support  the  larger  and  more  economically  and  recreational ly  important  ■■ 
species.  Macroinvertebrates  are  of  special  importance  in  stream  environments 
because  of  their  role  in  recycling  large  amounts  of  organic  detritus  introduced 
from  uplands. 

Macroinvertebrates  were  selected  as  a parameter  for  this  study  because 
they  are  sensitive  to  environmental  changes  and  thus,  are  important  indicators 
of  water  quality.  Natural  or  man-induced  fluctuations  in  the  physical - 
chemical  characteristics  of  a lentic  system  are  reflected  by  shifts  in  benthic 
conmunity  structure.  They  tend  to  remain  at  fixed  locations  and  have  a rela- 
tively short  life  span  of  usually  a year  or  less  and  therefore,  reflect  both 
the  present  and  recent  past  environmental  conditions. 

Macroinvertebrate  species  composition  (population  size,  number  of  taxa, 
and  diversity)  is  primarily  dependent  on  three  factors  --  water  quantity, 
water  quality,  and  substrate  composition. 

Water  quantity  limits  species  within  a site.  For  example,  some  organisms 
prefer  large,  deep  lakes  while  others  are  found  in  small,  shallow  lakes. 

Water  quality  is  also  a significant  factor  in  determining  the  assemblage  of 
macroinvertebrates.  Principal  parameters  include  oxygen,  temperature,  hardness, 
related  cations,  and  dissolved  solids.  The  most  important  of  these  is  oxygen. 
Increasing  evidence  shows  that  adverse  effects  can  be  expected  on  macroinverte- 
brates in  environments  supersaturated  with  oxygen  as  well  as  those  that  are 
undersaturated.  Aquatic  macroinvertebrates  are  also  affected  by  temperature 
extremes.  The  Aquatic  Life  Advisory  Committee  (1956)  indicates  that  benthic 
communities  in  temperature  zones  are  adapted  to  seasonal  fluctuations  of 
temperature  between  0 and  32°C  (32  - 90‘’F). 


Substrate  is  the  most  important  determinant  in  species  composition 
(Hynes,  1960).  There  is  a direct  relationship  between  the  amount  of  available 
surface  area  and  species  abundance  and  diversity.  That  is  to  say,  there  are 
more  hiding  and  foraging  places  in  a rock  or  pebble  bottom  than  in  a sand  or 
mud  bottom.  Contribution  of  organic  matter,  particularly  from  plants,  is  also 
an  important  consideration.  Aquatic  plants  increase  the  abundance  and  diver- 
sity of  benthic  organisms  viz.  there  is  more  surface  area,  periphytic  food 
organisms,  food  from  the  plants  themselves,  and  detritus  from  which  to  feed. 
Beck  (1954)  states,  "...after  careful  examination  of  many  streams,  diversity 
of  fauna  was  primarily  the  result  of  one  factor  --  the  diversity  of  habitat." 
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Environmental  Conditions 


The  winter  of  1976-1977  was  severe  and  unusually  low  temperatures  were 
recorded  throughout  the  season.  Ice  cover  on  Waconda  Bay  prevented  the 
commencement  of  Phase  II  until  late  February.  These  conditions  influenced 
the  composition  of  the  macroinvertebrate  community  by  reducing  populations. 

Low  lake  levels  caused  a significant  reduction  in  water  surface  area  in  the 
upper  end  of  Waconda  Bay.  Effluent  from  VAAP  was  confined  to  a relatively  narrow, 
shallow  channel  for  some  3,000  feet  downbay  after  reaching  reservoir  pool 
elevation.  Lack  of  mixing  in  this  channel  caused  significant  concentrations 
of  effluent  to  extend  further  downbay  than  would  be  expected  at  higher  pool 
elevations . 

Detritus  was  present  at  all  stations,  but  was  abundant  only  at  Station  6. 

A mat  of  filamentous  algae  grading  from  heavy  to  sparse  between  stations  3-6 
was  also  observed.  The  slender  spikerush  (Eleocharis  acicularis)  was  found  in 
small  quantities  at  Station  12. 

Table  14  presents  the  characteristics  of  the  bay  sediments  at  the 
various  stations.  Except  for  Stations  8,  9,  and  10,  all  samples  were  obtained 
at  comparable  depths.  During  Phase  II,  15  replicate  grabs  were  collected  at 
each  station.  In  addition,  14  Hester-Dendy  5-plate  units  were  placed  at  each 
station  location  for  colonization  over  a 4-week  period. 

Natural  Substrate  Colonization 


The  results  for  this  part  of  the  study  are  shown  in  Table  15  and  are 
expressed  as  numbers  of  organisms/square  meter.  The  various  taxa  are  grouped 
into  two  categories  representing  chironomids  and  other  forms.  Except  at  Station 
6, midges  composed  50  percent  or  more  of  the  total  population.  Three  forms 
were  most  abundant.  These  were  Procladius  sp,  Coelotanypus  spp.  and  Ch i ronoiiius 
sp.  Procladius  sp.  dominated  at  Stations  1-5  and  may  represent  forms 
more  tolerant  to  VAAP  waste. 

Total  population/m^  ranged  from  586  to  1,318.  Numbers  of  species  were 
generally  unifoniily  distributed  with  9 stations  showing  10  - 15  species  each. 

The  lowest  number  of  taxa  (7)  was  noted  at  Station  2, 

Other  forms  such  as  oligochaetes  were  not  abundant.  The  largest  population 
amounted  to  793  organisms/m^  at  Station  6.  Populations  of  fingernail  clams  were  rather 
sparse  between  Stations  1 and  4.  Further  downbay  where  effluent  effects  would 
be  less,  these  fingernail  clams  were  observed  ranging  from  60  to  190/m2.  Other 
organisms  demonstrated  similar  patterns.  Hexagenia  mayflies  were  found  at 
Station  3 indicating,  at  least,  some  tolerance  to  the  waste  in  this  location. 

A large  increase  in  these  insects  was  observed  at  Stations  4,  5,  and  6.  Shannon- 
Weaver  Diversity  is  expressed  both  for  culled  and  unculled  data  in  Table  16. 

For  natural  substrates  values  range  to  slightly  above  2 for  each  of  the  culling 
routines.  Station  3 represents  the  area  of  highest  diversity  followed  by 
Station  8.  Stations  1 and  2 in  the  maximum  zone  of  impact  were  the  least 
diverse.  In  general,  there  seemed  to  be  no  relationship  between  diversity  and 
station  location  vis-a-vis  the  waste  outfall. 
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CHARACTERISTICS  OF  THE  BENTHIC  SUBSTRATE  AT  THE  SAMPLING  STATIONS 
WITHIN  WACONDA  BAY  AND  THE  REFERENCE  BAY,  AS  OBSERVED  FROM 


VAAP  MACROINVERTEBRATES.  NATURAL  SUBSTRATE,  MARCH  1977 


Unculled  Data 


The  data  as  analyzed  by  cluster  analysis  are  presented  in  Figures  14  and 
15.  The  phenograms  are  representative  of  two  options  in  the  analysis  pro- 
cedure. In  Figure  14,  organisms  which  were  present  at  4 percent  or  less 
of  the  total  population  at  any  station  were  culled  and  no  longer  considered 
in  the  analysis  procedure.  Using  the  Pinkham-Pearson  Similarity  Index,  the 
option  of  mutual  absence  important  was  employed.  Figure  15  shows  the  clus- 
tering pattern  based  on  all  organisms  that  were  collected  and  where  mutual 
absence  was  considered  unimportant.  These  two  routines  were  run  to  examine  the 
impact  on  station  clustering  of  organisms  occurring  on  a seldom  or  rare  basis. 
As  the  figures  show,  the  clustering  is  relatively  similar.  In  both  cases 
stations  1 and  2 are  highly  unrelated  to  the  remaining  stations.  Stations  3, 

4,  and  5 are  somewhat  similar.  No  really  consistent  pattern  exists  for  the 
remaining  stations  indicating  that  stations  6 through  9 in  Waconda  Bay  are 
not  distinguishable  from  stations  10  through  12  in  the  reference  bay. 


Artificial  Substrate  Colonization 

Overall,  colonization  of  the  Hester-Dendy  units  was  sparse  and  densities 
were  greatly  reduced  over  those  in  the  benthos.  For  example,  chironomids  were 
2 to  3 orders  of  magnitude  less  than  those  populating  sediments  in  Waconda 
Bay.  Combined  data  from  five  Hester-Dendy  units  per  station  are  presented 
in  Table  1 7, 

Chironomids  dominated  artificial  substrates  at  Stations  7 through  9 
in  Waconda  Bay,  and  Stations  10  and  11  in  the  reference  bay.  Cricotopus 
spp.  accounted  for  about  two-thirds  of  the  Chironomidae.  Several  species  of 
Cricotopus  are  involved;  however,  presently  there  are  no  effective  means  to 
separate  these  species.  Four  species  of  Dicrotendipes  and  two  species  of 
Glyptotendipes  comprised  most  of  the  remaining  chironomids.  These  six  feed 
upon  planktonic  materials.  Cricotopus  was  present  throughout  Waconda  Bay,  and 
was  abundant  in  the  reference  bay.  Dicrotendipes  and  Glyptotendipes  were 
relatively  common  in  the  reference  bay  and  were  frequently  found  in  the  open 
bay  stations  (6  through  9)  of  Waconda  Bay.  Two  other  genera,  Ablabesmyia  and 
Tanytarsus , were  frequently  found  in  the  reference  bay,  but  rarely  found  in 
Waconda  Bay. 

Of  the  non-chironomid  fauna,  only  two  of  the  nine  found  on  artificial 
substrates  were  present  in  significant  numbers.  The  larvae  of  Caenis,  a 
small  herbivorous  mayfly,  reached  its  greatest  population  at  Station  6 where 
it  formed  61  percent  of  the  population.  Densities  fell  steeply  either  upbay 
or  downbay.  It  was  not  collected  at  Stations  1,  2,  9,  and  ID.  Dligochaetes 
occurred  sporadically  in  Waconda  Bay,  but  were  the  dominant  taxa  on  reference 
bay  artificial  substrates.  The  larvae  of  Argia  sp.,  a predaceous  damsel  fly, 
were  found  only  at  Stations  1 and  2.  This  genus  has  previously  been  reported 
only  from  unpolluted  situations  (Beck,  1954)  in  Florida,  but  Beck  (personal 
communication)  believes  this  may  not  hold  true  for  some  non-Florida  species. 

The  highest  densities  and  numbers  of  species  occurred  in  the  reference 
bay  stations  which  suggests  some  inhibition  in  all  Waconda  Bay  samples. 
Considering  total  number  of  species  per  station  the  stations  seem  to  group  as 
follows:  stations  1 to  5 averaged  3.3  species;  stations  6 to  9 averaged  9.5 
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Pinkham-Pearson  Biotic  Similarity 


FIGURE  14  . PHENOGRAM  OF  MACROINVERTEBRATES,  NATURAL 
SUBSTRATE,  CULLED  AT  4%,  MUTUAL  ABSENCE 
IMPORTANT.  COPHENETIC  CORRELATION 
COEFFICIENT,  0.917. 


Pinkham-Pearson  Biotic  Similarity 


FIGURE  15  . PHENOGRAM  OF  MACROINVERTEBRATES,  NATURAL 
SUBSTRATE,  UNCULLED,  MUTUAL  ABSENCE  UN- 
IMPORTANT. COPHENETIC  CORRELATION 
COEFFICIENT,  0.900. 


VAAP  MACRO IN VERTEBRATES,  ARTIFICIAL  SUBSTRATE,  UNCULLED,  MARCH  1977 


species;  stations  10  and  11  averaged  15.5  species.  Station  12  was  excluded 
because  the  sampling  rack  had  been  overturned  in  the  field  resulting  in  the 
Hester-Dendy  units  being  suspended  at  about  one  half  the  depth  of  the  units 
at  other  stations. 

Using  the  five  replicates  from  each  station,  an  analysis  of  variance 
was  performed  on  the  organism  density  data.  These  data  were  first  transformed 
to  In  (N  + 1).  The  calculated  F value  (degrees  of  Freedom  11,  48)  was  47.3 
which  is  significant  at  less  than  0.0001.  Using  Duncan's  multiple  range  test 
(5%  criteria)  to  compare  treatment  means  it  was  found  that  significant  differ- 
ences in  population  density  are  as  follows: 

1 , 2,  3 < 7,  5,  4,  9,  8 < 6 < 10  < 11  < 12 

Table  18  shows  diversity  values  for  the  community  colonizing  the  Hester- 
Dendy  units.  The  data  are  expressed  using  the  culled-unculled  format.  Treat- 
ment of  the  data  by  either  option  shows  Station  8 to  be  the  most  diverse  with 
a maximum  of  2.02.  The  lowest  diversities  were  observed  between  Stations  1 
through  5,  suggesting  that  these  stations  were  under  greater  stress. 

Cluster  analysis  using  the  Pinkham-Pearson  approach  was  performed  on 
data  from  the  Hester-Dendy  units  (Figures  16&17),  The  unculled  data  were  analyzed 
on  a mutual  absence  unimportant  basis  whereas  when  the  data  were  culled,  mutual 
absence  was  considered  important.  In  the  culled  data.  Stations  1 to  5 group 
separately  as  do  Stations  11  and  12.  In  the  unculled  data,  at  Stations  1,  3,  4 
and  5 are  grouped  together.  Station  2 is  unlike  any  other  station  and  Stations 
11  and  12  are  highly  similar.  Thus,  in  both  cases.  Stations  1 to  5 seem  to 
separate  from  the  remaining  group. 
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TABLE  18 

SHANNON  - WEAVER  DIVERSITY  INDICES  (H)  FOR 
VAAP  ARTIFICIAL  SUBSTRATES  MACROINVERTEBRATES 
MARCH,  1977 


Station 

Magnitude 

Culled 

Magnitude 

Unculled 

1 

0.69 

0.69 

2 

0.54 

0.54 

3 

0.95 

0.95 

4 

1.01 

1.01 

5 

0.36 

0.36 

6 

1.39 

1.57 

7 

1.55 

1.55 

8 

2.02 

2.02 

9 

1 .80 

1.80 

10 

1.69 

1.81 

11 

1.66 

1.80 

12 

1.11 

1.20 

Pinkham-Pearson  Biotic  Similarity 
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FIGURE  16  . PHENOGRAM  OF  MACROINVERTEBRATES,  ARTIFICIAL 
SUBSTRATE,  CULLED  AT  4%,  MUTUAL  ABSENCE 
IMPORTANT.  COPHENETIC  CORRELATION 
COEFFICIENT,  0.902 


Pinkham-Pearson  Biotic  Similarity 
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FIGURE  17  . PHENOGRAM  OF  MACROINVERTEBRATES,  ARTIFICIAL 
SUBSTRATE,  UNCULLED,  MUTUAL  ABSENCE  UN- 
IMPORTANT. COPHENETIC  CORRELATION 
COEFFICIENT,  0.882. 
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DISCUSSION  AND  CONCLUSION 


r 


In  relating  environmental  effects  to  munitions  concentrations  in  a 
field  study  such  as  this,  one  should  understand  the  several  assumptions  and/ 
or  limitations  that  are  inherent.  First,  the  concentrations  at  any  site 
vary  with  time.  For  the  natural  substrate  sampling,  the  concentrations 
measured  from  March  1-5  are  presumed  to  be  "normal"  or  "average"  and  the 
biologic  communities  are  assumed  to  be  in  equilibrium  with  these  concentra- 
tions. For  the  artificial  substrates,  measurements  were  made  of  munitions 
in  the  effluent  during  the  incubation  period  but  no  bay  samples  were  taken. 
However,  from  the  effluent  concentrations  and  an  understanding  of  how  the 
effluent  mixes  with  bay  waters,  one  can  predict  possible  resultant  bay 
concentrations.  Also,  the  plant  effluent  contains  a complex  mixture  of 
munitions- related  organic  compounds  plus  significant  concentrations  of  other 
ions,  particularly  nitrate  and  sulfate.  Thus  the  observed  effects  may  not 
be  due  to  only  the  five  specific  munitions  compounds  measured.  Even  with 
these  limitations,  certain  conclusions  can  be  relatively  firm.  Other 
conclusions  can  only  be  tentative  due  to  the  factors  noted  above  and 
because  the  effects  indicated  are  subtle. 

Overall,  the  data  appear  to  be  very  consistent.  The  water  chemistry, 
artificial  substrate  periphyton  population  density,  chlorophyll  a,  biomass 
and  viable  cell  ratios,  and  artificial  and  natural  substrate  macroinverte- 
brate community  structure  all  indicate  that  Stations  1-5  are  affected  by  the 
VAAP  discharge.  The  only  data  that  do  not  support  this  conclusion  are  sedi- 
ment chemistry  and  natural  substrate  periphyton.  Given  that  the  effluent  does 
appear  to  be  having  an  effect,  the  problem  thus  becomes  one  of  relating  the 
effects  to  specific  causes. 


First,  at  what  concentration  levels  were  no  measurable  effects 
observed?  In  the  natural  substrate  no-effect  coul  d be  discerned  in  either 
number  of  species,  density,  or  community  structure  for  both  benthic  and 
periphyton  organisms  at  stations  6 through  9 as  compared  to  reference  bay 
stations  10  through  12.  During  the  sampling  period  total  concentration  of 
the  five  munitions  compounds  averaged  less  than  25  gg/l  at  stations  6 
through  9.  In  the  December  sampling  average  munitions  levels  in  this  area 
ranged  as  high  as  27  ug/1.  Hence,  it  appears  that  for  total  munitions  con- 
centration less  than  25  yg/l  no  effects  would  be  expected. 

At  the  opposite  end  of  this  range,  effluent  munitions  concentrations 
during  the  incubation  period  for  the  artificial  substrate  samples  averaged 
about  575  pg/1  and  reached  a peak  of  2370  ug/1 • Concentrations  at  stations  1 
through  5 should  not  have  exceeded  these  values.  At  this  level,  definite 
effects  were  noted  in  periphyton  population  density  and  community  structure 
and  in  macroinvertebrate  numbers  of  species  and  community  structure.  It 
cannot  be  determined  whether  the  peak  or  average  concentration  was  most 
influential  in  the  result. 

Between  these  extremes,  effects  were  subtle  and  less  pronounced. 

In  the  five  day  period  prior  to  the  incubation  of  the  artificial  substrate, 
munitions  concentrations  in  the  upper  channel,  stations  1 through  5,  were 
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averaging  about  80  yg/l . Data  from  the  December  survey  and  results  of 
prior  studies  suggest  that  the  long  term  average  munitions  concentration 
in  this  area  might  actually  be  somewhat  higher.  Even  so,  the  only  observed 
effect  on  macro in vertebrate  infauna  was  a slight  shift  in  community 
structure.  Substrate  variation  in  periphyton  samples  may  have  masked  slight 
shifts  in  community  structure  in  that  compartment  also. 

During  the  incubation  of  the  artificial  substrate  samples,  it  is  pre- 
dicted that  the  concentration  of  munitions  compounds  at  stations  6-9  could 
have  gone  up  to  50  - 100  yg/l.  For  macroinvertebrates , density,  number  of 
species,  and  community  structure  at  stations  6-9  were  different  than  at 
stations  in  the  reference  bay.  For  periphyton,  the  numbers  of  species  seemed  to  be 
slightly  lower  at  stations  6 to  9 as  compared  to  the  reference  bay  stations. 

These  differences  were  not  as  substantial  as  those  observed  between  the  upbay 
and  reference  areas.  However,  they  were  of  such  magnitude  as  to  eliminate 
them  from  the  definite  no-effect  category. 

In  the  1975  study  (Sullivan,  et  aj^. , 1977),  biostimulation  was  observed 
in  certain  portions  of  Waconda  Bay.  In  this  winter  survey  no  such  biostimula- 
tion was  observed.  Possibly  the  low  temperature  encountered  prevented  this. 

In  summary,  the  work  reported  herein  indicates  no  environmental  effects 
on  periphyton  or  macroinvertebrates  from  a complex  TNT  manufacturing  effluent 
at  a total  munitions  concentration  of  less  than  25  yg/1.  Definite  effects 
were  noted  during  a period  when  total  munitions  averaged  500-600  yg/l  with 
a peak  of  over  2000  ug/1 . The  data  further  suggests  the  tentative  conclusion 
that  concentrations  in  the  range  of  50-100  pg/l  produce  minimal  effects.  Again 
it  should  be  noted  that  the  concentrations  of  other  non-munitions  unique 
compounds  changed  as  the  munitions  concentrations  changed.  These  changes 
probably  were  factors  also  in  the  biologic  changes  observed. 

These  results  are  in  good  agreement  with  the  conclusions  of  the  1975 
study  (Sullivan  et  1977).  In  that  work,  toxicity  was  indicated  in  the 
range  of  50-350  ug/l  total  munitions  and  no  effects  were  observed  at  concen- 
trations of  20  ug/1  and  less. 
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FIGURE  A-1.  RHODAMINE  B DYE  CONCENTRATION-SURFACE 

(ppb  by  weight) 


FIGURE  A-3.  CONDUCTIVITY-SURFACE  (ymhos/cm) 


FIGURE  A-4.  TOTAL  MUNITIONS  CONCENTRATION,  ug/1 . 


i’ 


APPENDIX  B 

ANALYTICAL  PROCEDURES 


i 


101 


I 


I 


LIST  OF  TABLES 

APPENDIX  B - ANALYTICAL  PROCEDURES 


table  DESCRIPTION 

{ B-1  SUMMARY  OF  ROUTINE  LABORATORY  ANALYTICAL  PROCEDURES 

FOR  WATER  SAMPLES 

B-2  SUMMARY  OF  ROUTINE  LABORATORY  ANALYTICAL  PROCEDURES 

FOR  SEDIMENT  SAMPLES 

( 


102 

^ — — _ 


ANALYTICAL  PROCEDURES 


Standard  Methods  or  EPA-approved  procedures  were  employed  to  characterize 
the  background  water  quality  (see  Table  B-1  for  specific  procedures).  Anmonia 
nitrogen  and  total  Kjeldahl  nitrogen  were  analyzed  on  a Technicon  Autoanalyzer 
II  using  the  manufacturer's  specified  methodology.  The  block  digestor  and 
the  sal icylate-nitroprusside  colorimetric  methods  are  currently  being 
researched  by  the  EPA  and  are  approved  for  use  in  NPDES  permit  monitoring 
(Gales,  EPA,  1977  personal  communication).  Analytical  procedures 
(Table  B-2  ) utilized  on  sediment  samples  were  primarily  from  Chemistry 
Laboratory  Manual  Bottom  Sediments  (EPA,  1969). 

Munitions  Analysis 

The  water  samples  were  collected  in  amber  glass  reagent  bottles  that 
were  pre-rinsed  in  acetone.  The  samples  were  refrigerated  until  analysis, 
which  consisted  of  extraction,  concentration  and  gas-liquid  chromatography. 
Sediment  samples  were  collected  in  glass  jars  that  had  been  pre-rinsed  with 
acetone.  Saran  wrap  was  used  as  a cap  liner. 

Extraction  of  Water  Samples,  December  1976.  A sample  of  250  ml  was 
measured  into  a clean  500  ml  separatory  funnel  equipped  with  a Teflon  stop- 
cock. Seventy-five  ml  of  ethyl  acetate  (pesticide  grade)  was  added,  the 
flask  stoppered,  and  shaken  for  2 to  3 minutes.  The  layers  were  allowed  to 
separate  and  the  lower  (water)  layer  drained  into  a second  500  ml  separatory 
funnel  and  again  extracted  with  50  ml  ethyl  acetate.  The  water  layer  was 
discarded.  The  extracts  were  combined  and  filtered  through  a plug  of  cotton 
previously  wetted  with  ethly  acetate.  The  separatory  funnels  were  rinsed 
with  an  additional  10  ml  of  ethyl  acetate  and  filtered  through  the  cotton 
plug.  The  ethyl  acetate  was  evaporated  to  a volume  of  2.5  ml  under  reduced 
pressure  with  the  flask  temperature  not  exceeding  40°C. 

Extraction  of  Water  Samples,  March  1977.  Same  procedure  except  500 
ml  samples  were  used. 

Extraction  of  Sediment  Samples.  In  order  to  dry  the  wet  sediments,  80 
gm  of  sodium  sulfate  was  added  to  20  gm  of  wet  sediment.  This  was  then 
packed  into  a chromatographic  column  and  extracted  for  one  hour  with  ethyl 
acetate.  The  extraction  was  followed  by  evaporation  of  the  ethyl  acetate 
extract  to  a volume  of  5.0  ml  under  the  same  conditions  as  described  earlier. 
Clean-up  techniques  are  discussed  in  the  recovery  studies  section  of  this 
procedure. 

Chromatography  of  Extracts,  December  1976.  Samples  were  chromatographed 
on  a 5 T72  ft.  x 1/8  in.  siTanized  glass  coTunin  packed  with  2%  OV  101/3"  QFl 
on  100/120  mesh  Gas  Chrom  Q.  A Varian  Model  1840  Gas  Chromatograph  with 
electron  capture  (EC)  and  Thermionic  (Alkali  Flame  Ionization  Detector) 

(AFID)  detectors  was  chosen.  The  readout  was  obtained  by  using  a Varian  Model 
285  Electronic  Integrator  which  was  recorded  permanently  by  a Beckman  1 mv, 

10  inch  scale  recorder.  Peak  areas  were  automatically  printed  by  integrator. 
Electron  capture  was  chosen  as  the  prime  detector  with  AFID  as  back-up  and 
confirmation  detector. 
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TABLE  B-1 


SUMMARY  OF  ROUTINE  LABORATORY  ANALYTICAL  PROCEDURES 
FOR  WATER  SAMPLES 


Parameter 


Procedure 


Chloride 
Total  Hardness 
Sul  fate 

Ammonia  Nitrogen 

Total  Kjeldahl  Nitrogen 


Nitrate-Nitrite  Nitrogen 


Standard  Methods . 408B;  Mercuric 
Nitrate  Method,  p.  304. 

Standard  Methods,  309B:  EDTA  Ti trimetric 
Method,  p.  202. 

Standard  Methods , 427C;  Turbidimetric 
Method,  BaCV^,  p.  496. 

Technicon  Corp.  Industrial  Method  #239- 
74  W/A:  Automated  colorimetric  salicy- 
late-nitroprusside  method,  1976. 

EPA,  STORET  #00625;  Acid  Digestion, 
Distillation,  Nesslerization , p.  175. 

or 

Technicon  Corp.  Industrial  Method  #376- 
75-W/A:  Acid  digestion  on  BD-20,  1975. 

Technicon  Corp.  Industrial  Method  #329- 
74-W/A:  Automated  colorimetric  salicy- 
late-nitroprusside  method,  1976. 

EPA,  STORET  #00630:  Automated  Cadmium 
Reduction  Method,  p.  207. 


EPA  1974,  Manual  of  Methods  for  Chemical  Analysis  of  Water  and  Wastes. 

Standard  Methods  for  the  Examination  of  Water  and  Wastewater  , Nth  Ed., 
1975,  APHA,  AWWA,  WPCF. 

Technicon  Industrial  Systems;  Tarrytown,  New  York. 
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TABLE  B-2 


SUMMARY  OF  ROUTINE  LABORATORY  ANALYTICAL  PROCEDURES 
FOR  SEDIMENT  SAMPLES 


Parameter 

Total  Kjeldahl  Nitrogen 
Nitrate-Nitrite  Nitrogen 

Total  Solids 

Total  Volatile  Solids 


Procedure 

Bottom  Sediments  - Great  Lakes;  Acid 
Digestion,  Distillation,  and  Titration 
with  0.02N  H^SO^,  p.  38. 

Bottom  Sediments  - Great  Lakes:  Acid 
Digestion,  p.  32. 

EPA,  STORET  #00630:  Automated  Cadmium 
Reduction  Method,  p.  207. 

Bottom  Sediments  - Great  Lakes: 

Gray  i met r i c Metliod , p.  85. 

BottOT  Sediments  - Great  Lakes: 
Gravimetric  Method,  p.  85. 


Chemistry  Laboratory  Manual  BottOT  Sediments,  EPA  1969,  compiled  by  Great 
Lakes  Region  Committee  on  Analytical  Methods. 

EPA  1974,  Manual  of  Methods  for  Chemical  Analysis  of  Water  and  Wastes. 


An  alternate  column  used  for  confirmatory  information  was  a 4 ft.  x 
1/7  in.  glass  column  packed  with  8 percent  UCW  98  on  80/100  mesh  Gas  Chrom 
Q.  Instrument  conditions  for  both  columns  and  detectors  were: 


Column  temperature: 
Injector  temperature: 
Detector  temperature: 
Carrier  gas: 
Electrometer  setting: 


l/B^C.  isothermally 

ZZO^C 

220‘'C 

Nitrogen  P 70  ml/min. 

10-10  afs  at  1 X attenuation  into 
integrator  with  appropriate  attenu- 
ation setting  for  recorder. 


Five  111  injections  of  sample  extracts  and  standards  were  first  injected 
onto  the  2%  OV  101/3%  AFl  column  using  the  EC  Detector.  Peaks  corresponding 
to  standards  were  noted  and  the  areas  compared.  Samples  and  standards  were 
next  injected  onto  USW-98  column  and  like  comparisons  were  made.  Likewise, 
samples  and  standards  were  injected  onto  the  2%  OV  101/3%  QFl  column  using  the 
Thermionic  or  Alkali  Flame  Ionization  Detector.  Again,  peaks  corresponding 
to  the  standards  were  noted  and  the  areas  were  compared.  Sample  peaks  which 
did  not  elute  at  the  same  times  as  the  standards  on  both  sets  of  columns  and 
detectors  were  rejected,  and  only  those  that  were  peaks  confirmed  on  both 
sets  were  quantitated.  Quantification  was  with  the  2%  OV  101/3%  QFl  column 
with  the  EC  detector. 


ChroMiatogr^hy  of  Extracts , Marc h^  1977.  Same  procedure  and  conditions 
as  above  with’ the  f'onowinV  exceptions;’ 

Electron  capture  detector  with  a standing  current  of  100%  at  16  x 10"^^ 

Attenuation  set  at  8 x 10-10 

Injection  volume  of  4.0  ul 

Column  temperature,  165"C 

Injector  temperature,  200°C 

Recovery  Studies,  DeceiPber  1976.  Recovery  studies  for  water  samples 
were  made  by  addition  of  stand’a’rds  to  tap  water  and  carrying  through  entire 
extraction,  concentration  and  GC  procedures.  Standards  were  prepared  in 
acetone  solution  prior  to  addition  to  water.  Previous  recovery  studies 
had  been  made  by  addition  of  standards  prepared  in  ethyl  acetate  solution, 
however,  it  was  thought  that  this  might  bias  the  study  since  ethyl  acetate 
would  be  used  for  subsequent  extraction.  Hence,  the  selection  of  water 
soluble  acetone  as  solvent  would  negate  any  positive  effect  which  the  standard 
solvent  miy  have  had. 

Since  previous  studies  had  shown  levels  of  less  than  1 ppm  TNT  and 
analogs  in  water  samples,  spiking  was  made  at  and  below  this  level  as 
follows : 

500  ml  water  was  used  in  all  cases. 
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COMPONENT 

QUANTITY  ADDED 

QUANTITY  RECOVERED 

% RECOVERY 

2,  4 DNT 

500  yg 

480 

96 

250  yg 

230 

92 

100  yg 

90 

90 

2,  6 DNT 

200  yg 

470 

94 

250  yg 

240 

96 

100  yg 

92 

92 

1,  3 DNB 

500  yg 

480 

96 

250  yg 

225 

90 

100  yg 

94 

94 

2,  4,  6 TNT 

500  yg 

460 

92 

250  yg 

230 

92 

100  yg 

95 

95 

SPIKING  STUDIES.  WATER  SAMPLES 


250  ml  sample  used 


SAMPLE 

TOTAL 

NO. 

COMPONENT 

yq  PRESENT 

yg  ADDED 

RECOVERY 

% RECOVERY 

AB  21 

2,4  DNT 

0 

25 

24 

96 

2,6  DNT 

0 

25 

21 

84 

TNT 

0 

25 

25 

100 

1,3  DNB 

0 

25 

22 

88 

AB  no 

2,4  DNT 

55 

25 

70 

88 

2,6  DNT 

24 

25 

45 

92 

TNT 

69 

25 

95 

101 

1,3  DNB 

1 

25 

22 

85 

AB  114 

2,4  DNT 

6 

25 

30 

97 

2,6  DNT 

3 

25 

25 

89 

TNT 

3 

25 

24 

86 

1,3  DNB 

0 

25 

23 

92 

AB  37 

2,4  DNT 

0 

62 

59 

95 

2,6  DNT 

0 

62 

63 

102 

TNT 

0 

62 

60 

97 

1,3  DNB 

0 

62 

58 

94 

AB  103 

2,4  DNT 

51 

62 

105 

93 

2,6  DNT 

26 

62 

81 

92 

TNT 

19 

62 

75 

93 

1,3  DNB 

1.5 

62 

60 

94 

AB  204 

2,4  DNT 

28 

62 

82 

91 

2,5  DNT 

15 

62 

73 

95 

TNT 

15 

62 

75 

97 

1,3  DNB 

1.5 

62 

58 

91 

107 


i 


w 


I 


Sediment  samples  were  column  extracted  with  ethyl  acetate  after  adding 
sufficient  anhydrous  sodium  sulfate  to  absorb  moisture  and  yield  apparently 
dry  product.  Samples  thus  extracted  contained  too  many  interfering  substances 
to  be  run  directly  on  GLC  and,  therefore,  were  subjected  to  florisil  column 
cleanup. 

Elution  of  florisil  with  ethyl  acetate  not  only  removed  the  TNT 
compounds,  but  most  of  the  interfering  substances  as  well.  It  was  found 
that  elution  of  florisil  column  with  6%  ethyl  ether  in  hexane  gave  quanti- 
tative recovery  of  the  compounds  of  interest  while  leaving  most  of  the 
extraneous  materials  on  the  column.  This  procedure  was  adopted  for  cleanup. 

Overall  recovery  averaging  approximate! v 80%  was  not  as  good  as  with 
water  samples.  A more  efficient  extraction  procedure  would  improve  recovery. 

Recovery  of  Standards  added  to  Anhydrous  Sodium  Sulfate 


COMPOUND 

QUANTITY  ADDED 

QUANTITY  RECOVERED 

% RECOVERY 

2,4  DNT 

100  ug 

94 

94 

2,6  DNT 

100  ug 

90 

90 

TNT 

100  pg 

95 

95 

2,4  DNT 

250  pg 

230 

92 

2,6  DNT 

250  pg 

125 

90 

TNT 

250  pg 

235 

94 

Indicated 

recovery  = 92%. 

Recovery  of  Standards  from 

Florisil  Cleanup  Procedure 

COMPOUND 

QUANTITY  ADDED 

QUANTITY  RECOVERED 

% RECOVERY 

2,4  DNT 

100  pg 

96 

96 

2,6  DNT 

100  pg 

95 

95 

TNT 

100  pg 

98 

98 

2,4  DNT 

250  pg 

240 

96 

2,6  DNT 

250  pg 

230 

92 

TNT 

250  pg 

240 

96 

Recovery  from  Spiked  Samples  of  Sediment 


SAMPLE 

NO. 

COMPONENT 

pg  PRESENT 

pg  ADDED 

TOTAL 

RECOVERY 

% RECOVERY 

W4-1 

2,4  DNT 

0 

25 

21 

84 

2,6  DNT 

0 

25 

20 

80 

TNT 

0 

25 

22 

88 

W4-10 

2,4  DNT 

3 

25 

22 

79 

2,6  DNT 

2 

25 

21 

78 

TNT 

5 

25 

26 

87 

W4-4 

2,4  DNT 

0 

50 

40 

80 

2,6  DNT 

0 

50 

38 

76 

TNT 

0 

50 

42 

84 

W4-6 

2,4  DNT 

0 

50 

43 

86 

2.6  DNT 

0 

50 

42 

84 

TNT 

3 

50 

42 

79 

Recovery  of  standards  added  to  samples  and  then  put  through  extraction 
and  cleanup  procedures  was  not  as  good  as  extraction  and  cleanup  without  the 
presence  of  sample,  thus  indicating  that  the  TNT  compounds  present  some 
difficulty  in  extraction  from  sediment  samples.  Overall  recovery  indicated 
recovery  is  approximately  80%. 

Recovery  Studies,  ^rch  1977.  Standards  were  prepared  in  acetone 
solution  by  diluting  stock  standard  (1  mg/ml)  with  acetone  to  100  iig/ml  or  0.1 
ug  per  pi.  The  combination  standard  then  contained  O.lng/pl  of  each  of  the 
following: 

1,  3 DNB 

1 , 3 , 5 TNB 

2,  4 ONT 
2,  6 DNT 

2,  4,  6 TNT 

Five  hundred  ml  of  distilled  water  was  added  to  each  of  three  separatory 
funnels  and  combination  standard  was  added. 

100  1 (10  pg)  to  No.  1 

500  1 (50  pg)  to  No.  2 

2.5  ml  (250  pg)  to  No.  3 

Contents  of  each  funnel  was  mixed,  the  funnel  stoppered  and  allowed  to  stand 
1 hour.  Contents  of  each  funnel  was  then  extracted  and  taken  through  same 
cleanup  concentration  and  chromatography  as  were  samples. 


Recoveries  were: 


Component 

Quantity 

Added 

Quantity 

Recovered 

. . . pg  . 

Percent  Recovered 

1.3  ONB 

10 

9.1 

91 

50 

46.8 

94 

250 

240 

96 

1,3,5  TNB 

10 

9.0 

90 

50 

47.0 

94 

250 

235.0 

94 

2.4  DNT 

10 

9.5 

95 

50 

43.0 

86 

250 

230.0 

92 

2,6  DNT 

10 

9.3 

93 

50 

48.0 

96 

250 

232.0 

93 

2,4,6  TNT 

10 

9.2 

92 

50 

51.1 

102 

250 

242.0 

97 

Studies  were  also  made  with  spiked  samples  similar  to  above  except 
that  known  increments  of  standards  were  added  to  actual  samples.  Two 
hundred  ml  samples  were  spiked  with  varying  amounts  of  standards  using  same 
procedures  as  with  distilled  water.  These  were  extracted,  concentrated, 
run  through  cleanup  and  chromatographed  same  as  original  samples  with  the 
following  results: 


Component 

Amount  originally 

1,3  DNB 

1,3,5  TNB 

2,4  DNT 

2,6  DNT 

2,4,6 

found  (ug) 

3.1 

12.5 

- 

1.7 

14.1 

Amount  added  (iig) 

10.0 

10.0 

10.0 

10.0 

10.0 

Total  present  (ug) 

13.1 

12.5 

10.0 

11.7 

24.1 

Assay  (ug) 

12.3 

12.3 

9.1 

10.5 

22.2 

Recovery 

93.9 

98.4 

91.0 

89.7 

92.1 

Component 

Amount  originally 

1,3  DNB 

1 ,3,5  TNB 

2,4  DNT 

2,6  DNT 

2,4,6 

found  (ug) 

- 

8.8 

1.1 

2.2 

22.6 

Amount  added  (ug) 

10.0 

10.0 

10.0 

10.0 

10.0 

Total  present  (ug) 

10.0 

18.8 

11.1 

12.2 

32.6 

Assay  (ug) 

9.3 

18.0 

10.3 

11.1 

31.2 

T.  Recovery 

93.0 

95.7 

92.8 

91.0 

95.7 

Component 

An)ount  originally 

1,3  ONB 

1,3,5  TNB 

2,4  DNT 

2,6  DNT 

2,4,6 

found  (ug) 

- 

0.2 

- 

6.5 

1.0 

Amount  added  (ug) 

5.0 

5.0 

5.0 

5.0 

5.0 

Total  present  (ug) 

5.0 

5.2 

5.0 

11.5 

6.0 

Assay  (ug) 

5.1 

4.8 

4.7 

10.7 

5.8 

% Recovery 

102.0 

92.3 

94.0 

93.0 

96.7 

no 


Component 

1,3  DNB 

1,3,5  TNB 

3,4  DNT 

2,6  DNT 

2,4,6  TNT 

Amount  originally 
found  (ng) 

. 

. 

1.5 

0.7 

1 .5 

Amount  added  (pg) 

5.0 

5.0 

5.0 

5.0 

5.0 

Total  present  (pg) 

5.0 

5.0 

6.5 

5.7 

6.5 

Assay  (pg) 

4.9 

4.8 

6.2 

5.6 

6.6 

% Recovery 

98.0 

96.0 

95.4 

98.2 

101.5 

Component 

1 ,3  ONB 

1,3,5  TNB 

3,4  DNT 

2,6  DNT 

2,4,6  TNT 

Amount  originally 
found  (pg) 

_ 

. 

20.8 

15.0 

5.9 

Amount  added  (pg) 

10.0 

10.0 

10.0 

10.0 

10.0 

Total  present  (pg) 

10.0 

10.0 

30.8 

25.0 

15.9 

Assay  (pg) 

9.7 

9.5 

27.9 

23.2 

15.1 

% Recovery 

97.0 

95.0 

90.6 

92.8 

95.0 

Sediment  samples  were  column  extracted  with  ethyl  acetate  after  adding 
sufficient  anhydrous  sodium  sulfate  to  absorb  moisture  and  yield  an 
apparently  dry  product.  Samples  thus  extracted  contained  too  many  interfering 
substances  to  be  run  directly  on  GLC  and,  therefore,  were  subjected  to 
florisil  cleanup,  eluting  with  6%  ether  in  hexane.  The  areas  of  the  plate 
containing  the  compounds  of  interest  were  then  removed,  extracted  with  ethyl 
acetate.  This  cleaned  up  sample  was  then  adjusted  in  volume  and  quantitated 
using  GLC. 


Results  of  spiking  studies  are  given  below; 


Component 

1,3  DNB 

1 ,3,5  TNB 

2,4  DNT 

2,6  DNT 

2,4,6  TNT 

Amount  originally 

present  (pg) 

1.1 

6.4 

- 

0.8 

4.7 

Amount  added  (pg) 

2.5 

2.5 

2.5 

2.5 

2.5 

Total  present  (pg) 

3.6 

8.9 

2.5 

3.3 

7.2 

Assay  (pg) 

3.0 

7.0 

2.0 

2.5 

6.0 

% Recovery 

83.3 

78.7 

80.0 

75.8 

83.3 

Component 

1,3,5  TNB 

2,4,6  TNT 

Amount  originally 

present  (pg) 

15.5 

8.6 

Amount  added  (pg) 

5.0 

5.0 

Total  present  (pg) 

20.5 

13.6 

Assay  (pg) 

16.5 

11.0 

% Recovery 

80.5 

80.9 

Component 

1,3,5  TNB 

2,4,6  TNT 

Amount  originally 

present  (pg) 

6.5 

6.2 

Amount  added  (pg) 

5.0 

5.0 

Total  present  (pg) 

11.5 

11.2 

Assay  (pg) 

9.0 

9.5 

% Recovery 

78.3 

84.8 

Because  of  several  factors  such  as  extraction  efficiency,  additional 


I cleanup  steps  and  handling  losses,  overall  recovery  from  sediment  samples 

i averages  about  80^  whereas  recovery  from  water  samples  was  considerably 

better. 
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D-1  PERIPHYTOMETER  SUSPENSION  UNIT  123 


METHODOLOGY  FOR  PERIPHYTON  COLLECTION 


Within  the  study  area  periphyton  samples  were  collected  on  natural 
and  artificial  substrates.  Forms  colonizing  natural  substrates  were  collected 
at  the  end  of  February,  while  those  on  exposed  glass  slides  were  not  collected 
until  the  end  of  March.  For  this  reason,  the  data  are  not  comparable  except 
in  a general  qualitative  sense. 

Three  samples  representative  of  natural  substrates  were  obtained  at 
each  site  and  preserved  in  5 percent  formalin.  These  were  oxidized  with 
hydrogen  peroxide  and  potassium  di chromate  to  remove  organic  matter.  The 
treated  sample  was  dried  on  an  18  x 18  coverslip,  mounted  in  Hyrax  (refrac- 
tive index  1.7),  and  examined  under  oil  immersion  (lOOOX).  For  qualitative 
comparisons,  a minimum  of  500  diatom  valves  were  identified  per  sample.  Raw- 
count  data  were  transferred  to  a coding  sheet.  This  information  was  key 
punched  prior  to  computer  analysis.  Several  routines  were  utilized  including 
culling  the  data  to  eliminate  species  which  occurred  only  rarely.  The 
collection  of  periphyton  on  glass  slides  is  considered  to  give  more  reliable 
station-to-station  comparisons  than  natural  substrates  because  of  influence 
by  substrate  specificity.  Standard  1 x 3 inch  microscope  glass  slides  were 
placed  in  periphyton  samplers  (Periphytomer  II)’^  one  inch  below  the  water 
surface  at  nine  stations  located  in  Waconda  Bay,  and  three  stations  located 
in  the  reference  bay  (Figure  3).  Peri  phytometers  were  placed  in  suspension 
units  as  shown  in  Figure  D-1 . Samplers  were  placed  according  to  recoiiinended 
procedures  outlined  in  Standard  Methods  for  the  Examination  of  Water  and 
Wastewater  14th  Edition,  ~(APHA,  T975T,  aTiTBioloqicaT  Field  and  Laboratory 
Methods  for  Measuring  the  Quality  of  Surface  Water  and  Effluents,  (Weber, 
1573). 


Peri  phytometers  were  anchored  in  a northerly  direction  utilizing  the 
suspension  unit  designed  for  this  study,  A total  of  8 racks  (64  slides) 
were  placed  in  the  suspension  units  at  each  station.  As  part  of  a supple- 
mental statistical  study  to  measure  variability  between  replicate  slides 
and  to  determine  the  best  practical  method  for  processing  periphyton  slides, 
additional  6-rack  (48  slides)  suspension  units  were  placed  at  Stations  1,  3, 

8,  9,  11  and  12.  Anchorline  adjustments  were  made  for  anticipated  lake  level 
fluctuations  of  2 feet  or  less.  All  slides  were  collected  over  a 2-day  period 
from  March  31  - April  1,  1977,  following  an  incubation  period  of  30  days. 
Sampling  devices  were  recovered  at  all  stations  except  9 and  10,  where  peri- 
phytometer  units  had  been  overturned  exposing  slides  to  the  air.  However, 
at  station  9,  the  statistical  6-rack  unit  of  slides  was  still  intact  and  was 
used  as  a substitute  for  the  vandalized  8-rack  units. 

At  the  end  of  the  30-day  incubation  period,  slides  from  both  the  8-rack 
and  6-rack  units  were  recovered  utilizing  the  following  collection  procedures. 
End  slides  (Positions  1 and  8)  from  the  8-rack  units  were  removed,  dried,  and 
placed  into  labeled  glass  bottles.  Slides  for  biomass,  chlorophyll  a,  diatom, 
filamentous  algae,  and  vital  stain  counts  were  collected  from  slide  positions 
2-7  using  the  random  number  generator  of  a Texas  Instruments  SR-51A 
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calculator.  The  priority  for  selecting  slides  from  each  rack  was: 


2 end  slides  - Air-dried 

2 slides  - Biomass  determination 

2 slides  - Chlorophyll  a 

1 slide  - Diatom  counts 

]_  slide  - Filamentous  algae  or  vital  stain  counts 

8 slides  - Total 


For  the  8-rack  periphyton  suspension  units  the  minimum  numbers  of  slides  for 
analyses  were: 


Analyses 


Minimum  Number  of  Replicates 
per  Station 


Chlorophyll  a^  16 
Biomass  16 
Diatom  counts  5 
Filamentous  algae  2 
Vital  stain  counts  ^ 

Total  42 


Slides  from  the  6-rack  suspension  units  were  all  air-dried  for  diatom 

counts. 


Chlorophyll  a Determination 

Except  for  Station  10,  a total  of  16  replicate  periphyton  slides  were 
collected  at  each  station  for  chlorophyll  a analysis.  While  in  the  field, 
periphyton  slides  were  placed  in  50  ml  of  a 90  percent  acetone  (v/v),  10  per- 
cent of  a saturated  MgCO^  solution  and  immediately  stored  in  the  dark  on  dry 
ice.  Prior  to  analysis, "^chlorophyll  was  extracted  for  24  hours  in  the  dark 
at  4"C.  To  facilitate  extraction,  slides  were  scraped  and  the  acetone  suspen- 
sion ground  30  seconds  at  500  rpm  in  a Potter-type  tissue  homogenizer. 

Following  extraction  chlorophyll  a^,  corrected  for  phaeophytin,  was 
determined  fluorometrical ly  after  the  methods  of  Yentsch  and  Menzel  (1963), 
Holm-Hansen,  ^ (1965),  Lorenzen  (1967),  and  Moss  (1968),  using  a Turner 

Design  Model  10  fluorometer.  Fluorometric  determination  of  chlorophyll  depends 
on  red  fluorescence  emitted  by  the  chlorophyll  a^ molecule  when  excited  by 
ultraviolet  light  and  is  100  times  more  sensitive  than  spectrophotometric 
analysis.  The  method  is  limited  to  chlorophyll  a^  only,  chlorophyll  ^ and  £ 
cannot  be  determined. 

The  chlorophyll  £ reference  solution  was  a purified  spinach  chlorophyll 
standard  (Product  No.  C5753,  Sigma  Chemicals,  St.  Louis,  MO)  calibrated  by 
spectrophotometric  chlorophyll  analysis. 

Acidification  of  chlorophyll  £ converts  it  quantitatively  to  phaeophytin, 
Reading  the  fluorescence  before  and  after  adding  one  drop  of  IN  HCl  to  the  sam- 
ple cuvette  allows  calculation  of  an  acid  factor  related  to  the  interference. 
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Periphytic  chlorophyll  ^ was  calculated  as  follows: 

Chlorophyll  a , ,(F)(r)(Ca)(inl  extract) 

(r-1  )[substrate  area  (mm)  10”'^] 

where:  Ca  = fluorometer  reading  before-fluorometer  reading  after  acidification 

_ standard  before  acidification 
standard  after  acidification 

p _ r Ca 1 r dilution  ration  fluorometer  -■ 

^fluorometer  reading-*  *-dilution  ratio  spectrophotometer-* 

Organic  Biomass. 

Sixteen  replicate  biomass  slides  were  collected  from  each  station.  In 
the  laboratory  each  slide  was  rehydrated  for  15  minutes,  accumulated  material 
scraped  from  the  slide  into  a graduated  cylinder,  and  then  resuspended  in  a 
total  volume  of  50  ml  distilled  water.  An  aliquot  of  the  suspension  was 
filtered  on  a tared,  fi red-glass  filter  (Gelman,  6FA),  the  ash-free  dry  weight 
determined  (APHA,  1975)  and  converted  to  grams  of  organic  matter  per  square 
meter. 

Diatom  Cell  Densities. 

At  least  five  replicate  slides  per  station  were  examined  for  diatom 
community  structure  and  cell  density  estimates.  Periphyton  growth  was  scraped 
from  glass  slides  into  tall  labeled  beakers  using  razor  blades  and  a rubber 
policeman.  Samples  were  oxidized  with  20  ml  of  50  percent  hydrogen  peroxide 
and  approximately  50  mg  of  potassium  dichromate.  The  solution  was  cooled, 
allowed  to  settle  for  24  hours,  decanted,  and  brought  to  a volume  of  50  ml. 
Preliminary  examination  of  oxidized  material  from  Stations  6-12  indicated 
the  samples  were  too  concentrated  to  count  accurately.  Therefore,  4-  and  8- 
fold  dilution  was  necessary  to  permit  observations  of  10  - 15  organisms  per 
microscope  field.  In  contrast,  samples  at  Stations  1 - 5 were  too  sparse  and 
had  to  be  concentrated  by  a factor  of  5 to  permit  accurate  counting. 

Permanent  slides  were  prepared  by  pipetting  0.4  ml  of  the  "cleaned" 
material  onto  an  18  x 18  coverslip  (324  mm2)  and  allowing  the  sample  to  dry 
at  65°C  (150°F)  on  a laboratory  hot  plate.  The  dried  coverslip  was  placed 
on  a standard  microscope  slide  containing  one  drop  of  Hyrax  mounting  medium 
(refractive  index  1.7)  and  the  slide  gently  heated  to  drive  off  the  toluene 
solvent.  When  cooled,  the  permanent  slide  was  labeled  with  station  number, 
date,  location,  and  dilution  factor.  Under  an  oil  immersion  lens  (Zeiss 
microscope,  lOOOX)  diatoms  were  identified  and  enumerated  to  the  species 
level,  where  possible,  utilizing  the  following  standard  taxonomic  references: 
Hustedt,  1930,  1962;  Schmidt,  et  al.,  1974-1959;  Huber-Pestalozzi , and  F. 
Hustedt,  1949;  and  Patrick  and  ReTiner,  1966,  1975.  Voucher  diatom  slides  were 
sent  to  Dr.  Charles  Reimer,  Academy  of  Natural  Sciences  of  Philadelphia,  for 
taxonomic  verification. 

Cell  densities  were  estimated  by  performing  field  counts  at  6-slide 
coordinates  randomly  selected  on  the  coverslip.  A total  of  10  microscope 
fields  were  examined  at  each  of  the  6 selected  coordinates  for  a total  of  60 
fields.  At  lOOOX,  each  microscope  field  represented  an  area  of  0.0182  iTin2 
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X 60  fields).  There- 


2 2 

with  a total  area  examined  of  1.092  mm  (i.e.  0.0182  mm 
fore,  cell  densities  were  determined  as  follows: 


Cells/mm  = 


Total  Area 
of  Coverslip 

Diatom  » (324  mm^) 

Counts  * Total  Area 
Examined 
(1.092  mm2) 


Original  Volume 
of  Periphyton 
Suspension 
(50  ml) 

Volume  of  Sample 
dried  on  coverslip 
(0.4  ml) 


Dilution 
Factor  / 


Original 
Surface  Area 
of  SI ide 
(3,871  mm2) 


Filamentous  Algae.  Two  replicate  periphyton  slides  per  station  were 
collected  from  artificial  substrate  units  located  in  Waconda  Bay  and  the 
reference  bay.  Slides  were  mechanically  scraped  with  razor  blades  and  pre- 
served in  5 percent  formalin.  Clumped  periphyton  material  was  dispersed  by 
blending  samples  in  a Waring  blender.  Algal  filaments  were  identified  using 
a 50  ml  plankton  sedimentation  chamber  and  a Zeiss  invertoscope  0.  Identifi- 
cations were  carried  to  species  level,  where  possible,  utilizing  the  following 
standard  references:  Drouet  (1968);  Prescott  (1962);  and  Desikachary  (1956). 

Vital  Stain  Counts.  Recent  studies  have  shown  that  natural  and  artificial 
substrate  periphyton  populations  contain  significant  numbers  of  diatom  frustules 
that  are  either  empty  or  have  less-than-robust  chloroplasts  (Pryfogle  and 
Lowe,  1976).  Standard  diatom  cleaning  procedures  using  ammonium  persulfate, 
acid  or  hydrogen  peroxide  oxidation  do  not  distinguish  between  diatoms  having 
metabol ical ly  active  chloroplasts  and  those  having  empty  frustules.  Under 
certain  conditions  an  investigator  may  wish  to  determine  the  metabolic  state 
of  diatoms  in  response  to  a pollution  source.  The  following  experimental  vital 
stain  procedure  was  conducted  at  11  stations  located  in  Waconda  Bay  and  the 
reference  bay  using  tetrazolium  violet.  Tetrazolium  salts  have  been  used  to 
determine  the  dehydrogenase  activities  of  bacterial  populations  in  activated 
sludge  and  to  a limited  extent  in  the  analysis  of  plankton  populations 
(Armitage,  1977).  The  enzymatic  reduction  of  the  tetrazolium  salt  (yellow) 
to  triphenyl  formazin  precipitate  (violet  stain)  visually  tags  diatom  cells 
with  active  cytochrome  systems.  As  a result,  cell  counts  can  be  qualitatively 
differentiated  into  one  of  three  metabolic  categories; 


Metabolic  Activity 


Cell  Description 


Viable  or  active 
Senescent 

Dead  or  fossil  forms 


Violet  precipitate  present  throughout 
cell,  or  concentrated  in  mitochondria. 
Chlorophyll  present,  no  violet  precipi- 
tate formed. 

No  violet  precipitate  or  chlorophyll 
present  in  frustule. 


Three  replicate  periphyton  slides  per  station  were  removed  from  the 
samplers  and  placed  in  glass  bottles  (75  ml)  containing  lake  water  collected 
at  the  site.  These  samples  were  then  placed  in  the  dark  and  1 ml  of  a 0.2 
percent  (w/w)  solution  of  paraiodotetrazol ium  was  added  to  the  sample  bottle. 
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The  slides  were  incubated  in  the  light  for  a period  of  12  hours  after  which  they 
were  preserved  with  3 percent  formalin.  Vital  stain  counts  were  made  utilizing 
a 50  ml  plankton  sedimentation  chamber  with  a Zeiss  inverted  microscope  equipped 
with  an  oil  immersion  lOOOX  lens.  Diatom  populations  were  differentiated  into 
metabolically  active,  senescent,  or  dead  (fossil). 


i 

I 
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STATISTICAL  METHODS 


I.  Introduction 

The  statistical  analyses  of  the  artificial  substrate  periphyton  data 
considered  the  following  conditions; 

- comparison  of  variability  in  density  and  diversity  estimates 
between  fields  on  a single  slide  with  variability  between  aliquots. 

- comparison  of  variability  in  density  and  diversity  estimates 
between  fields  and  aliquots  with  variability  between  slides. 

- comparison  of  variability  in  density  and  diversity  estimates  when 
one  slide  is  used  with  that  when  a composite  of  five  slides  is 
used. 

- comparison  of  the  eleven  sampled  stations  with  respect  to  diversity 
and  density. 

- at  one  station  a procedure  for  determining  dominant  species  to  a 
specified  degree  of  precision  was  determined. 

Shannon's  measure  of  diversity 
s , 

H = -I  In  Pi, 
i = l 

' th  - 

was  used  where  p^  is  the  proportion  of  the  i species  in  the  sample  (p^  = 

s is  the  number  of  observed  species,  and  we  use  natural  logarithms. 

Basharin  (1959)  has  shown  that  H is  a biased  estimate  of  the  population  diversity 

S 

H = p.  In  p. 
i = l ^ ^ 

where  p^  is  the  true  proportion  of  species  i in  the  population,, and  S is  the 
true  number  of  species  present  in  the  population.  The  bias  of  H is  negative, 
implying  that  R tends  to  underestimate  H.  Basharin  shows  that 


Bias 


s-1 

2N 


where  N is  the  total  number  of  organisms  observed.  Thus  a bias-adjusted 
estimate  of  H 

Hfi  = ^ 

can  be  obtained.  Basharin  also  derives  the  standard  error  of  Hg  showing  that 

5 

S.E.(Hg)  = Pi  - 
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which  can  be  estimated  by 

s:E.(Hg)  = 1n2  p.  - Hg2] 

This  allows  one  to  measure  the  reliability  of  the  diversity  estimate  and  to 
compare  the  station  diversities  statistically. 

Turning  to  density,  the  logarithmic  transformation 

L.  = ln(N.  + 1) 

was  used  to  transform  N.,  the  total  number  of  organisms  per  observational 
unit  (usually  a group  of  fields)  to  a variable  L-  which  more  nearly  satisfies 
the  assumptions  necessary  for  an  analysis  of  variance.  Specifically,  the 
variable  will  often  have  a distribution  which  is  approximately  normal  when 
is  skewed  to  the  right.  This  transformed  density  measure  can  now  be  used 
to  make  statistical  comparisons  of  the  various  sampling  procedures,  and  of 
the  stations. 


II.  Comparison  of  Diversity  and  Density  Variability  Between  Fields  With 
That  Between  Aliquots 

A.  Description  of  Experiment 

One  slide  was  selected  from  each  of  four  stations,  3,  8,  11,  and 
12.  Ten  aliquots  were  taken  from  each  slide,  and  120  fields  or  500 
organisms,  whichever  came  first  (with  a 60  field  minimum),  were 
counted  on  the  microscope  slide  made  from  each  aliquot.  Counts  were 
made  in  units  of  10  fields.  This  experiment  enables  one  to  compare 
the  variability  between  fields  to  that  between  aliquots  at  each  station. 

B.  Diversity 

As  shown  in  Section  VIII-A,  one  may  analyze  the  bias-adjusted 
diversities  using  a nested  analysis  of  variance  model.  Table  E-1 
presents  the  results  of  the  analyses  for  the  four  stations. 

C.  Density 

The  same  nested  analysis  of  variance  model  which  was  used  for 
diversity  was  used  to  analyze  L,  the  logarithmic  transform  of 
density.  The  results  are  shown  in  Table  E-2. 

D.  Conclusion 

Whether  the  replicates  consist  of  groups  of  fields  from  one 
microscope  slide,  or  of  fields  from  different  aliquots,  one  can 
pool  the  results  when  estimating  diversity  and  density.  There  is 
no  evidence  that  the  aliquot  variability  and  the  field  within  aliquot 
variability  differ. 
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TABLE  E-1 


COMPARISON  OF  DIVERSITY  VARIABILITY  BETWEEN  ALIQUOTS 


TO  THAT 

BETWEEN  FI 

[ELDS  FOR  4 

STATIONS 

Station 

Source 

DF 

ss 

MS 

F 

3 

A1 iquot 

9 

3.4132 

.3792 

1.00 

F ields(Al iquot ) 

110 

41.826 

. 3802 

H 

A1 iquot 

9 

0.1648 

.01831 

0.46 

Fields(Al iquot ) 

60 

2.4041 

.04007 

11 

Aliquot 

9 

0.1731 

.01923 

0.86 

Fields{Al iquot) 

50 

1.1240 

.02248 

12 

A1 iquot 

9 

0.3569 

. 03966 

1.06 

Fields(Al iquot) 

70 

2.6068 

.03724 

None  of  the  F values  are  sign! ficant  at  even  the  .10  level.  There  is  no 
evidence  that  aliquoting  introduces  variability  into  the  estimate  of 
diversity. 
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TABLE  E-2 

COMPARISON  OF  DENSITY  VARIABILITY  BETWEEN  ALIQUOTS 
TO  THAT  BETWEEN  FIELDS  FOR  4 STATIONS 


I 

i 


ation 

Source 

OF 

Si 

MS_ 

F 

3 

Aliquot 

9 

1.9445 

.2161 

.81 

Fields(Aliquot) 

110 

29.3589 

.2669 

8 

A1 iquot 

9 

1.1033 

.1226 

.70 

Fields(Aliquot) 

60 

10.5679 

.1761 

11 

A1 iquot 

9 

1.4106 

.1567 

1.01 

Fields(Al iquot) 

50 

7.7954 

.1559 

12 

A1 iquot 

9 

1.7829 

.1981 

.98 

Fields(Al iquot) 

70 

14.1604 

.2023 

As  with  diversity,  there  is  no  evidence  to  suggest  that  aliguoting  contributes 
variability  to  the  estimate  of  density. 


III.  Comparison  of  Diversity  and  Density  Variability  Between  Slides  To 
That  Within  Slides 

A.  Experiment 

Nine  slides  were  selected  from  each  of  the  stations  3,  8,  11,  and 
12,  and  a count  was  obtained  for  60  fields  (10  fields  at  each  of  6 
random  locations)  on  a sinijl®  microscope  slide  prepared  from  each 
experimental  slide.  Thus  the  variability  between  slides  can  be 
determined,  and  then  compared  to  that  within  slides. 

B.  Diversity 

The  variance  estimate  of  diversity  between  slides  is  based  on  the 
60  field  counts  of  the  single  microscope  slide;  while  the  within 
slide  variability  is  estimated  from  the  10  field  counts  from  the 
previous  experiment.  Since  the  same  stations  are  sampled,  one  would 
expect  approximately  six  times  more  organisms  in  60  fields  than  in 
10,  which  theoretically  increases  the  diversity  variance  by  a factor 
of  six.  An  appropriate  test  statistic  to  compare  these  variance 
components  is  therefore 

c _ 6(Variance  Between  Slides) 

(Tariance  i/Tt)iin  ^lidesl 

The  results  of  this  test  are  summarized  in  Table  E-3. 


The  five  percent  level  of  F is  about  two  in  each  case.  It 
appears  that  the  between  slide  variability  is  larger  than  the  within 
slide  variability  at  least  at  stations  B and  11.  That  is,  there  is 
evidence  that  the  diversity  variance  between  slides  may  be  slightly 
larger  than  that  within  slides. 


C.  Density 

The  assumption  that  the  logarithmic  transform  L of  density  is 
approximately  normal  plays  a key  role  in  the  comparison  of  density 
variability  between  and  within  slides.  Using  standard  results  for 
the  lognormal  distribution,  one  can  find  an  estimate  of  within  slide 
variance  to  be 


2vl  + 
e 


-1) 


2 

where  y.  is  the  mean  of  the  log  transform  L and  S . Js  the  v.iriance 
of  L.  ^The  same  formula  can  be  used  with  the  between  slide  data  to 
calculate  S^n,  the  estimate  of  between  slide  variability.  Then  one 
can  test  that  the  between  slide  variability  exceeds  the  within  slide 
va  i a b i 1 i ty  by  compu  t i ng 


r 
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where  the  factor  of  6 is  present  because  6 times  more  fields  are 
involved  in  each  density  count  for  the  between  slide  calculations. 

The  results  of  the  analysis  are  given  in  Table  E-4.  Note  that  this 
F test  is  only  approximate,  since  the  distributions  are  log-normal 
rather  than  normal . 

D.  Conclusion 

There  appears  to  be  slightly  more  variability  between  slides 
than  within  slides  for  both  diversity  and  density.  The  implication 
is  that  one  should  obtain  replication  by  sampling  several  slides 
rather  than  aliquoting  a single  slide.  This  will  protect  against 
the  possibility  that  the  sampled  station  has  more  variability  between 
than  within  slides,  as  some  of  the  experimental  ones  did.  The  number 
of  slides  which  must  be  sampled  will  be  considered  in  Section  VI. 

IV.  Comparison  of  Diversity  and  Density  Variability  When  One  Slide  is 
Used  With  That  When  Five  Slides  Are  Composited 

A.  Experiment 

Five  slides  were  composited  and  one  microscope  slide  formed,  on 
which  60  fields  were  counted.  Five  replicates  of  this  experiment  were 
performed  at  each  of  four  stations,  3,  8,  11  and  12. 

B.  Diversity 

In  order  to  determine  whether  the  compositing  reduces  variability, 
one  can  compare  the  composited  variance  to  the  single  slide  variance 
from  the  previous  experiment  (see  III-A).  The  results  are  summarized 
in  Table  E-5. 

C.  Density 

The  same  comparison  can  now  be  made  on  the  logarithmic  transforms 
of  the  densities  (see  Table  E-6)- 

D.  Conclusion 

The  compositing  of  five  slides  does  not  seem  to  reduce  the  variability 
of  either  diversity  or  density  estimates.  Since  more  replication  (and 
thus  better  variance  estimates)  can  be  obtained  by  analyzing  the  slides 
separately,  the  single  slide  method  is  probably  preferable  to  compositing 
slides. 


V.  Station  Comparisons 

A.  Experiment 

Since  the  previous  results  indicated  that  a replicated,  single 
slide  experiment  at  each  station  was  preferred,  this  experiment  was 
performed  at  seven  other  stations,  five  replicates  at  each  station. 

B.  Diversity 

The  stations,  their  bias-adjusted  diversity,  and  the  standard 
error  (see  Basharin,  1959)  are  given  in  Table  E-7. 
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TABLE  E-4 

COMPARISON  OF  DENSITY  VARIABILITY  BETWEEN 
SLIDES  TO  THAT  WITHIN  SLIDES 


ation 

Source 

DF 

Variance 

F 

3 

Between  Slides 

8 

250.05 

2.99 

Within  Slides 

119 

13.95 

8 

Between  Slides 

8 

17,489.1 

3.42 

Within  Slides 

69 

853.06 

11 

Between  Slides 

8 

9,293.4 

1.10 

Within  Slides 

59 

1,409.7 

12 

Between  Slides 

8 

12,394.2 

2.35 

Within  Slides 

79 

880.88 

There  is  ample  statistical  evidence  to  conclude  that  the  density  variability 
between  slides  exceeds  that  within  slides  at  stations  3,  8,  and  12. 


TABLE  E-5 


COMPARISON  OF  COMPOSITED 

DIVERSITY 

VARIANCE 

TO  SINGLE  SLIDE 

VARIANCE 

Station 

Source 

DF 

Variance 

F 

3 

Single 

8 

.07152 

1.21 

Composite 

4 

.05920 

8 

Single 

8 

.01731 

1.46 

Composite 

4 

.01185 

11 

Single 

8 

.01068 

3.29 

Composite 

4 

.00325 

12 

Single 

8 

.009864 

0.84 

Compos i te 

4 

.01170 

Since  the 

F value  for  significance  at  the 

0.10  level 

is  F = 3.94 

, no 

evidence 

that  composite  diversity  variance 

1 exceeds  single  slide 

variance 

is  found. 


I 


i 


1 
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TABLE  E-6 


COMPARISON  OF  COMPOSITED  DENSITY  VARIANCE 
TO  SINGLE  SLIDE  VARIANCE 


Station 

Source 

DF 

Variance 

F 

3 

Single 

8 

.1129 

15.6 

Composite 

4 

.007224 

8 

Single 

8 

.1236 

2.85 

Compos i te 

4 

.04330 

11 

Single 

8 

.03815 

0.38 

Composite 

4 

.1017 

12 

Single 

8 

.04572 

0.62 

Composite 

4 

.07419 

No  consistent  pattern  in  these  variance  comparisons  is  found.  There  seaiis 
to  be  significantly  more  variance  in  the  single  slide  experiment  than  the 
composite  slide  experiment  at  Station  3,  but  somewhat  less  at  Station  11. 
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TABLE  E-8 

POPULATION  DENSITIES  FOR  ALL  STATIONS 

Density  (number/mm^) 

37,293 
27,803 
22,462 
19,698 
18,936 
14,911 
750 
174 
144 
85 


Thus  the  conclusion  is  that  the  station  densities  are  aligned  as  follows: 

11  > 9,  12,  7,  6,  5,  4,  1,  3,  2 
8 > 6,  5,  4,  1,  3,  2 
9,  12,  7,  6 > 5 > 4 > 1,  3 > 2 


! 


Station 

11 

8 

9 

12 

7 

6 

5 

4 

1 

3 

2 
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TABLE  E-7 

DIVERSITIES  FOR  ALL  STATIONS 


Station  Diversity  Standard  Error 


6 

0.716 

0.035 

7 

1.189 

0.033 

5 

1.608 

0.043 

8 

2.045 

0.025 

9 

2.175 

0.032 

4 

2.440 

0.074 

11 

2.557 

0.020 

1 

2.573 

0.074 

12 

2.790 

0.019 

2 

2.918 

0.089 

3 

2.929 

0.061 
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An  analysis  of  variance  comparison  was  cleaned  inappropriate  due  to 
the  probable  inequality  of  the  diversity  variances.  Thus  the  stations 
were  compared  in  a pairwise  fashion.  The  lines  on  the  left  of  the 
station  numbers  reflect  the  results  of  comparing  the  diversities  by 
using  the  two  sample  Z statistic,  connecting  stations  which  cannot  be 
declared  significantly  different.  The  formula  for  Z is 

1 - (Station  i Diversity  - Station  j Diversity) 

>^(Standard  Error  i)^  + (Standard  Error  j)^ 

Since  Basharin  (1959)  shows  that  Shannon's  diversity  estimate  is  asymp- 
totically normal,  this  statistic  should  approximate  a standard  normal 
random  variable  when  the  true  population  diversities  are  the  same,  and  it 
does  not  necessitate  equal  variances.  '<Jhen  the  value  of  Z is  large 
(negative  or  positive),  the  implication  is  that  the  true  station  diver- 
sities differ.  A conservative  .01  significance  level,  Z = 2.58,  was  used 
since  so  many  comparisons  were  made.  Table  E-7  shows  that  the  stations' 
true  diversities  are  arranged  as  follows: 

6 < 7 < 5 < 8 < 9 < 4,  11,  1 < 12,  2,  3 
C.  Density 

The  logarithmic  density  transforms  were  subjected  to  a standard  analysis 
of  variance.  The  station  comparison  yielded  a highly  significant  F value 
of  319.81,  leaving  no  doubt  that  the  station  densities  differ.  The 
transformed  densities  were  then  subjected  to  Duncan's  (1975)  multiple  com- 
parison procedure  to  determine  which  stations  differ.  The  estimated 
densities  (untransformed)  are  shown  in  Table  E-8,  with  the  solid  lines 
on  the  left  of  the  station  numbers  connecting  those  stations  which  one 
cannot  conclude  as  different  at  the  .05  significance  level. 

VI.  Sample  Size  Determination 

Perhaps  the  most  crucial  decision  which  must  be  made  in  planning  an  experi- 
ment is  the  sample  size,  or  number  of  replicates  to  be  processed.  One  can  use 
information  gleaned  from  previous  experiments  to  plan  more  carefully. 

In  Section  III  it  was  concluded  that  replication  of  slides  was  more  desirable 
than  replication  of  aliquots  from  a single  slide.  Hence,  it  was  decided  that  the 
effect  on  a single  slide  consist  of  counting  ten  adjacent  fields  at  each  of  six 
randomly  selected  locations  (60  total  fields).  The  question  then  becomes:  how  many 
slides  should  be  sampled? 

A.  Diversity 

Since  diversity  is  analyzed  on  an  untransformed  scale,  it  is  easier 
to  work  in  absolute  rather  than  percentage  differences.  The  following 
equation  (Steel  and  Torrie  1960)  is  for  the  absolute  difference  between 
two  stations'  diversities  which  can  be  detected  with  90  percent  confidence 
for  n replicate  slides: 

D = 2(1.645)a^  = 4.653i^^/  vlT 
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TABLE  E-8 

POPULATION  DENSITIES  FOR  ALL  STATIONS 


Station 

11 

8 

9 

12 

7 

6 

5 

4 

1 

3 

2 


Densi ty  ( number /nw^ ) 

37,293 

27,803 

22,462 

19,698 

18,936 

14,911 

750 

174 

144 

85 

28 


Thus  the  conclusion  is  that  the  station  densities  are  aligned  as  follows: 

11  > 9,  12,  7,  6,  5,  4,  1,  3,  2 
8 > 6,  5,  4,  1,  3,  2 
9,  12,  7,  6 > 5 > 4 - 1,  3 > 2 
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where  oi  is  the  standard  deviation  of  the  diversity  based  on  single 
slide  estimates.  This  assumes  that  the  stations  being  compared  have 
the  same  standard  deviation  oj.  The  probability  of  falsely  concluding 
a difference  exists  is  set  at  a = .10. 

A plot  of  D versus  n for  the  four  experimental  stations  is  shown 
in  Figure  E-1.  Note  that  for  all  but  Station  3,  where  the  density  of 
organisms  was  very  low,  one  can  detect  a diversity  difference  of 
slightly  less  than  0,3  with  90%  confidence  using  five  replicate  slides. 


B.  Density 

Since  the  statistical  analysis  of  density  is  performed  on  a 
logarithmic  scale,  the  differences  one  can  detect  between  station 
densities  will  be  expressed  in  percentages.  The  absolute  percentage 
difference  which  can  be  detected  with  90%  confidence,  with  a = .10, 
is  given  by 

P = 100{exp[4.653oj/  - 1} 

where  0,  is  the  standard  deviation  of  the  transformed  density  for 
single  slide  estimates.  A plot  of  p versus  n for  the  four  experimental 
stations  is  shown  in  Figure  E-2.  Note  that  one  must  have  at  least  five 
replicate  slides  before  he  can  be  90%  sure  of  detecting  density 
differences  of  100%.  It  is  obvious  that  density  comparisons  require 
more  replication  than  do  diversity  comparisons  to  achieve  similar 
rel lability. 


C.  Conclusions 


In  order  to  determine  the  necessary  sample  size  for  estimating  and 
comparing  station  diversities  or  station  densities,  a preliminary  estimate 
of  the  variance,  a,  , must  be  obtained.  One  procedure  is  to  count  approxi- 
mately 60  fields  on  a single  slide  from  the  station  of  interest,  and  then 
to  estimate  by  the  following. 

Diversity: 


1 


s 


p,  In  p. 


' 2 

- VI 


(See  Section  I for  notation  definitions) 


Then  the  number  of  slides  necessary  to  detect  a difference,  D,  in 
diversities  with  90  percent  confidence  for  two  stations  having  approximately 
the  same  standard  deviation  is  given  by 

n = 21.65  d^^/D^ 


If  the  stations  of  interest  have  very  different  variances,  compute  the 
average  variance  to  use  in  the  formula  for  n. 

Density:  The  computation  of  an  estimated  variance  for  the  density  is 
more  difficult.  Perhaps  the  best  method  is  to  sample  3 slides  from  each 
station,  count  60  fields  on  each  slide,  and  then  calculate  the  sample 
variance,  s,'^,  for  the  3 log-transformed  density  observations.  Then  the 
total  number  of  slides  which  need  to  be  sampled  in  order  to  detect  a 


of  Replicate  Slides, 


FiajPEE-2-  SW'PLE  SIZE  TO  «EE’  OE’iSITY  CRI^EPIA 


TABLE  E-9 

DOMINANT  SPECIES  FOR  STATION  8,  PERIPHYTON 
ARTIFICIAL  SUBSTRATE 


percentage  difference  P between  two  stations  is  given  by 

n = 21.65  s^^/[ln(l+P/100)l^ 

2 

Again,  s,  may  be  obtained  by  averaging  the  various  station  variances 
if  preliminary  samples  are  obtained  for  several  stations. 

VII.  Determination  of  Dominant  Species 

An  analysis  of  dominant  species  was  conducted  using  9 replicate  slides 
at  Station  8.  The  basis  for  dominant  species  selection  was  that  each  species 
make  up  at  least  3%  of  the  pooled  sample  at  the  station.  Six  species 
satisfied  this  criteria  at  Station  8.  The  data  are  shown  in  Table  E-9. 

The  first  question  to  be  answered  was:  Can  one  pool  all  nine  re- 
plicates and  treat  each  species  composition  estimate  as  a binomial  estimate. 
For  example,  the  total  count  of  species  1 in  all  replicates  was  258,  which 
is  7.9%  of  the  total  count  of  3264.  If  one  can  use  binomial  theory,  the  95T. 
confidence  interval  for  species  1 percent  composition  is 

7,90  i 2 X /IMMim 


or 


7.90  + 0.94 

Thus,  one  would  estimate  with  95?.  confidence  that  between  6.96  and  8.84 
percent  of  the  population  consists  of  species  1. 

A contingency  table  test  was  conducted  to  deteniiine  whether  there 
is  significant  difference  between  replicates.  The  table  consisted  of  7 rows 
corresponding  to  the  6 dominant  species  and  all  others  combined  and  9 columns 
corresponding  to  the  9 replicate  slides.  The  X'  statistic  of  110.8  is  signi- 
ficant at  the  .0001  level  when  compared  to  the  distribution  with  43  degrees 
of  freedom.  Thus,  it  is  concluded  that  a significant  difference  exists 
between  the  replicate  estimates  of  percent  composition  for  the  6 species. 

Since  the  replicates  are  all  taken  from  the  same  population  (station),  the 
percent  composition  of  each  replicate  estimates  the  true  composition  at  the 
station.  The  fact  that  the  x*  test  shows  a significant  difference  in  repli- 
cate percentages  warns  that  the  pooling  and  use  of  binomial  theory  is  inadvisable. 

Using  the  theory  given  in  Section  VIII-B,  we  can  obtain  a sample 
estimate  of  variability  without  assuming  the  binomial  distribution  applies. 

Both  the  binomial  and  sample  estimates  are  shown  in  Table  E-10. 

It  is  obvious  that  if  one  were  to  use  the  binomial  formula  for  variance, 
the  result  would  be  an  underestimate.  Thus,  in  determining  the  sample  size 
necessary  to  estimate  the  percent  composition  for  a species  with  a specified 
accuracy,  the  binomial  formula  will  underestimate  the  necessary  sample  size. 

Several  methods  can  be  employed  to  more  realistically  determine  the 
sample  size  necessary.  For  example,  suppose  one  wanted  to  estimate  the 
percent  composition  attributable  to  species  3.  Diatpma  t(^nue  v.  eJjpnjjatum 
(from  Table  E-9),  the  most  dominant,  to  within  2%  with  95?  confidence. 

Although  the  binomial  bound  is  already  less  than  2?  (1.72")  with  3,264 
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TABLE  E-10 

VARIANCE  AND  BOUND  BASED  ON  BINOMIAL  THEORY  AND  ACTUAL  SAMPLING  RESULTS 


Species  Estimate  Variance  95%  Confidence  Bound 


Binomial 

Sample 

Binomial 

Samp!  1 

1 

7.90 

.0727 

.1203 

0.94 

1.21 

2 

3.28 

.0317 

.0401 

0.62 

0.70 

3 

40.35 

.2407 

1.0856 

1.72 

3.64 

4 

18.60 

.1514 

.7033 

1.36 

2.94 

5 

3.83 

.0368 

.0509 

0.67 

0.79 

6 

12.93 

.1126 

.3205 

1.17 

1.93 

organisms  counted,  the  more  realistic  sample  estimate  is  3.64^,  The 
number  of  organisms  necessary  to  reduce  this  bound  to  Zt  is  found  by 
solving; 


P _ 4 (Estimated  Variance)  ^ iq4 
( Bound 

= X 10^  = 10.856 

(2)^ 

This  works  out  to  about  21  more  replicate  slides  at  the  present  rate  of 
approximately  363  organisms  per  slide  (60  fields). 

The  main  disadvantage  of  the  above  method  is  that  one  must  have  an 
estimated  sample  variance  before  n can  be  calculated.  Thus,  several  replicates 
(at  least  three  would  be  advisable)  must  be  analyzed  before  the  necessary  sample 
size  can  be  estimated.  A cruder  alternative  is  available  when  preliminary 
sampling  cannot  be  conducted.  First,  one  can  estimate  the  sample  size  using 
the  binomial  formula,  then  multiply  by  a "safety  factor"  to  adjust  for  replicate 
variability.  The  appropriate  factor  is  the  ratio  of  the  true  variance  to  the 
sample  variance.  In  our  example,  TableE-lft  the  factor  ranges  from  about  1.5 
for  the  species  at  the  lower  end  of  the  dominant  scale  to  about  4 for  the 
species  with  high  percent  composition. 

For  example,  if  one  were  to  use  the  binomial  variance  to  estimate 
the  sample  size  necessary  to  estimate  the  species  3 percent  composition  to 
within  2%,  with  95%  confidence,  one  would  get: 

n a 4 (Binomial  Variance)  iq4 
(Bound)^ 

^ X 10^  = 2407 

(2)2 

Using  the  safety  factor  of  4,  since  species  3 is  at  the  upper  end  of  the 
dominant  percentages,  one  finds  that  about  9600  organisms  should  be  counted. 

This  is  not  too  far  from  the  estimate  of  10,856  obtained  using  the  sample 
results. 

In  summary,  if  slide  replicates  are  to  be  used  for  estimating  percent 
compositions,  the  variability  of  the  estimates  may  be  considerably  more  than 
is  expected  assuming  a binomial  distribution.  More  realistic  estimates  of 
variability  can  be  calculated  from  the  actual  sample  results,  and  these  in 
turn  can  be  used  to  estimate  sample  sizes  necessary  fora  specified  accuracy 
of  estimation.  If  the  binomial  distribution  is  used  to  determine  sample  size, 
a safety  factor  should  be  used  to  adjust  for  replicate  variability.  To  be 
safe,  this  factor  should  be  at  least  2 for  species  at  lower  percentages 
(5-10%),  and  3 to  4 for  species  above  10%. 
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VIII.  Detailed  Procedures 

A.  Analysis  of  variance  for  Shannon's  Diversity 

Suppose  one  is  estimating  a station's  true  diversity  by  selecting  a 
single  slide,  preparing  n.  aliquets  (and  thus  n^^  microscope  slides),  and 
then  counting  n^  sets  of  TO  fields  each  on  each  slide  (see  section  II  of  this 
report  for  the  analysis  of  such  an  experiment).  A statistical  model  for  such 
an  experiment  would  be: 


“jk  ' “ * "j  * ‘jk 


where 


th 


is  the  Shannon  diversity  estimate  for  the  k group 
of  10  fields  in  aliquot  j 
H is  the  true  Shannon  diversity, 

A.  is  a random  effect  associated  with  aliquot  j 

is  a random  effect  associated  with  the  k group  of 
10  fields  on  aliquot  j. 


Assume: 


where 


E(Aj)  = 0,  Var(Aj)  = o^^/N^. 

E(e Jk).  0,  »ar(Ej,). 

Cov  (A.,  . 0 

number  of  organisms  in  k^''  group  of  10  fields 
in  aliquot  j 
he 

l(?j  of  organisms  in  aliquot  j 


Thus 


E(H.k)  = H 

Var  t oe^/Nj^ 


Basharin  (1959)  shows  that, 

Var  (H..  ) = [I  Piln^  Pi  - /N. 


jk 


so  that  this  model  is  in  agreement  to  the  order  of  Njj^. 
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Now  one  can  apply  the  usual  analysis  of  variance  methods  to  this  nested 


model  for  First  calculate 


A *^c 

H - .1,  .1.  H.,/n,n^ 


j = 1.  2, 


j=l  k=l 

In  terms  of  the  model. 


so  that 


H.  = H 


H = H 


* ^ I I 


Var(Hj) 


n 2 2 n^ 

O , n_  C 


y i 

^3  Y k^-1  'Jk 


1 ^-l 


Where  N.  = ( J tj~) 

•J  k=l  '^jk 


2 'ifl 

o*  A 


— 1 
Var(H)  = — jj-  I ij—  + 


"a  "c 


1 


— 7 ^ JT  — ? — 7 ^ ^ N 

V j V^c  J'' 


"a1 


“7~2  ■ 

'’A  "c  =N 


where 


N = ( I . N = (I  f N 

j=l  ■ j k ^-k 


1 )-l 
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Define 


SSA 


" H)2 


= n I (Hj  - H) 
'j=l  J 


Thus 


= I [(H.  - H)  - (H  - H)]^ 
^ j=l  J 


n_  I (H.  - H)^  - n n (H  - H)‘ 
c J c « 


EfSSA)  = -£^  * ^ 
2 2 

"c°A 


= (n,-l)[-^  - -f-l 


EtMSA}  = Et^SSA}  = ^ 


2 ^ 2 
e 

A-  ="A"c 


Secondly, 


'’a  "c 


SSC  = 


77  x2 


I I (H.,^  - H ) 
j=l  k=l  J 


"a  "c 


= I I t(Hjk  - H)  - (H,  - H)]‘ 


j=l  k=l 


"A  "c 


= .1.  .1.  (^^jk  ■ - "r  .1.  (H.  - H)‘ 


j=l  k=l 


j = l 


/n_-lvo 

E(SSC)  = -4^ 


E(MSC)  = £{7;^  "It,  ■ SSC}  = 


V? 


153 


If  one  forms  the  approximate*  F ratio, 

c MSA 
^ ■ MSC 

2 

he  gets  a test  of  the  hypothesis  Hoi'^A  = 0,  or  that  the  aliquoting  contributes 
no  extra  variability  (beyond  that  of  the  sets  of  fields  within  aliquots)  to  the 
diversity  estimate.  Large  values  of  F would  imply  that  > 0,  or  that  the 
aliquoting  does  contribute  additional  variability.  These  tests  are  conducted 
in  section  II.  B. 


Note  that  this  same  reasoning  could  be  extended  to  additional 
hierarchies.  For  example,  one  could  similarly  analyze  an  experiment  involving 
stations,  slides  within  stations,  and  fields  within  slides  within  stations. 

In  this  way  one  can  determine  which  variance  components  seem  to  contribute 
most  to  the  variability  of  the  diversity  estimate.  These  tests  were  conducted 
in  sections  II  to  IV  of  this  report. 

B.  Procedure  for  Estimation  of  Percentage  of  Specific  Organism 

Assume  r replicate  slides,  with  n,,  np,  n organisms 
identified  on  each  slide.  Let  y,,  yo,  ....  y„  represent  the  number 
of  a particular  species  in  each^repiicate.  If  p is  the  true  fraction 
of  the  population  which  is  the  species  of  interest,  then  one  may 
estimate  p by  the  r fractions 


Yl  . y2  - Vr 

Pi  ■ nj’  P2  ■ 02*  ^r  “ n^ 


Central  limit  theory  implies  that  these  sample  Tractions  will 
be  approximately  normally  distributed,  with  mean  p.  If  the  binomial 
assumption  is  made,  the  variance  of  p.  is  p(l~P) . In  the  periphyton 

n. 

sampling  experiment,  there  is  evidence  that  this  value  underestimates 
the  actual  variance,  hence  it  is  proposed  that  the  assumption 

2 

Variance  (p. ) = ^ be  used 

where  is  an  unknown  constant,  exceeding  p(l-p).  One  can  estimate 
by  the  following  sample  statistic: 


I n^(p^-p)^ 

i=l  ^ 




where  p is  the  pooled  estimate  of  p,  i.e. 


r r 

I y.  I y 

- = i=l  ^i 

r R 

I n 
i=l  "i 


♦Basharin  (1959)  proves  that  the  diversity  estimate  is  asympotically  normal, 
which  lends  credence  to  the  assumption  that  this  ratio  possesses  an  F distribution. 
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Then  use  p to  estimate  p,  but  instead  of  estimating  the  variance  of 
p by  p(l-p)  , use 

^ ' s^ 

Var(p)  = ^ 


To  estimate  a percentage  composition  to  within  a specified  bound 
B and  with  a confidence  level  of  C percent,  first  find  the  z value 
from  the  normal  tables  corresponding  to  C,  i.e.  z = 2 for  C = 95  percent, 
z = 1.65  for  C = 90  percent,  etc.  Then  solve 


z2s2 

n = -|- 


X 10" 


4 

where  the  10  factor  is  used  of  B is  given  as  a percent,  but  eliminated 
if  B is  given  as  a decimal  fraction. 


IX.  Biologic  Interpretation  of  Results 

The  statistical  analysis  has  shown  that  a significant  source  of 
variability  in  determining  either  diversity  or  density  on  artificial 
substrate  periphyton  samples  is  the  difference  between  fields  on  an 
individual  slide.  In  the  procedure  utilized  here,  the  slides  were 
prepared  by  placing  a 0.4  ml  aliquot  on  an  18  x 18  nm  coverslip.  A 
random  number  procedure  was  used  to  determine  the  location  on  the  cover- 
slip  to  conduct  the  actual  count.  Ten  adjacent  fields  were  counted 
at  each  selected  coordinate,  i.e.  if  60  total  fields  were  counted, 
six  locations  on  the  coverslip  were  utilized. 

The  data  for  the  separate  10  field  units  shows  a considerable 
variation  in  total  numbers.  Locations  near  the  edge  of  the  coverslip 
tend  to  have  much  lower  organism  density  than  locations  near  the  center. 
The  reason  for  this  is  that  when  the  0.4  ml  aliquot  dries,  it  does  so 
from  the  edge  to  the  center.  This  has  the  effect  of  concentrating 
organisms  at  the  center  of  the  coverslip.  Hence,  there  is  a very  uneven 
distribution  of  diatoms  that  is  particularly  apparent  when  a random 
number  procedure  is  utilized  to  determine  the  counting  position  on  the 
coversl ip. 

Given  the  variation  described  above,  it  is  possible  that  the 
significance  of  aliquot  to  aliquot  or  slide  to  slide  variation  may  have 
been  masked.  Hence,  if  the  variation  due  to  drying  of  the  aliquot  on  the 
coverslip  could  be  reduced  then  the  aliquot  to  aliquot  and  slide  to  slide 
variation  could  be  assessed  in  a more  precise  manner.  A re-examination 
of  these  variables  would  be  appropriate. 

An  obvious  needed  improvement  in  processing  technique  is  the 
development  of  a procedure  that  does  not  concentrate  organisms  near  the 
coverslip  center.  Lacking  this,  it  is  important  that  a randomized  proce- 
dure be  used  to  select  coverslip  locations  which  counteracts  the  natural 
tendency  to  count  near  the  center  where  most  of  the  organisms  are  con- 
centrated. 
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APPENDIX  F 

MACRO INVERTEBRATE  METHODS 


LIST  OF  FIGURES 


METHODOLOGY  FOR  MACROINVERTEBRATE  COLLECTION 


Aquatic  macroinvertebrates  were  collected  from  natural  and  artificial 
substrates  in  March,  1977,  at  12  stations.  Natural  substrates  were  sampled 
with  a petite  Ponar  dredge.  Hester-Dendy  artificial  substrates  were  sus- 
pended approximately  1.5  - 3.0  feet  below  the  surface.  Fifteen  replicates 
of  the  natural  substrate,  and  fourteen  units  of  artificial  substrates,  were 
collected  at  each  station  to  minimize  natural  variability. 

To  decrease  variability  a special  unit  was  fabricated  for  Hester- 
Dendy  plate  suspension  (Figure  F-1).  This  sampling  device  was  designed  and 
built  by  Mr.  B.C.  Pruitt  of  Water  and  Air  Research,  Inc.  The  unit  consists 
of  PVC  pipe  molded  into  a circle  with  a circumference  of  approximately  2-1/2 
feet.  Fourteen  5 plate  units  are  attached  to  the  suspension  frame  by  wing 
nuts.  Floats  are  attached  to  a line  of  predetermined  length.  The  entire 
apparatus  is  anchored  at  selected  station  locations.  For  purpose  of  this 
survey, units  were  suspended  18  inches  below  the  surface. 

In  the  field,  dredge  samples  of  the  natural  substrate  were  washed  in  a 
bucket  sieve  (U.S.  Standard  No.  30  mesh)  and  bottled.  Rose  Bengal  dye  was 
then  added  to  facilitate  laboratory  sorting.  Samples  were  preserved  in  10 
percent  formalin.  Natural  substrate  samples  were  rewashed  in  the  laboratory 
and  picked  in  a white  enamel  pan  partially  filled  with  water.  After  sorting, 
organisms  were  placed  in  vials  containing  95  percent  ethanol.  Chironomid 
larvae  were  mounted  in  polyvinyl -lactophenol  for  microscopic  identifications. 

Samples  from  all  stations  were  analyzed  and  identifications  were  made 
to  the  lowest  practical  taxonomic  level.  Fifteen  natural  substrate  replicates 
were  processed  for  stations  1,  5,  and  11,  and  five  were  processed  for  the  re- 
maining nine  stations.  Five  of  the  14  artificial  substrate  replicates  (whole 
units)  were  processed  for  all  twelve  stations.  The  community  structure  indices 
which  were  computed  include  the  Shannon-Weaver  diversity  index  and  the  Pinkham- 
Pearson  index  of  biotic  similarity. 

Organism  Loss  From  Processing 

A problem  encountered  with  the  quantification  of  benthic  macroinver- 
tebrates was  the  loss  of  organisms  during  laboratory  processing.  This  occurs 
when  benthic  samples  are  rinsed  in  a 30-mesh  sieve  to  remove  formalin,  dye, 
and  the  finer  sediments.  Therefore,  we  saved  the  rinse  water  and  examined 
it  for  macroinvertebrates  for  the  five  replicates  at  Station  1 in  order  to 
estimate  losses  from  this  technique.  This  experiment  was  unnecessary  for 
artificial  substrates  since  they  were  not  rinsed  in  a sieve. 
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The  table  below,  with  the  five  replicates  pooled,  enumerates  the  results: 


Taxa 


Number  of  Organisms  Retained  By: 
30-mesh  Sieve  Rinse  Pan 


Procladius  sp. 

Tanypus  neopunctipennis 

Oligochaeta 

Nema  toda 


88 

7 

1.5 

1 


3 

3 

0.5 

0 


Unidentified  Coleopteran 


0 


Totals 


98.5  6.5  >:=105 

(93. 8X)  (6.2%) 


A similar  experiment  was  conducted  on  a leaf  pack  sample.  Of  an 
estimated  400-600  organisms  within  the  sample,  approximately  100-200  passed 
through  the  sieve.  Almost  all  those  lost  were  chironomids  and  we  feel  that 
a 6.2  percent  loss  is  acceptable.  The  loss  is  apparently  only  partially 
related  to  organism  size  when  considering  those  with  small  diameters. 
Experience  has  shown  that  there  is  very  little  loss  when  samples  contain 
moderate  to  large  amounts  of  filamentous  algae  or  sand.  Apparently,  sand 
plugs  many  of  the  sieve  openings  and  filamentous  algae  entangle  the  organisms 
which  results  in  higher  sieve  retention. 


The  sediment  in  the  study  area  was  composed  almost  entirely  of  clay 
and  silt.  This  matrix  was  of  a pudding-like  consistency.  It  was  typically 
orange  or  brown  (oxidized)  on  top,  becoming  gray  (reduced)  underneath.  The 
top  orange  layer,  as  observed  from  dredged  benthic  grabs,  was  about  1-2  cm 
at  Stations  1 and  2.  It  was  2-4  cm,  at  Stations  3,  4 and  5,  but  reduced  to 
a thin  veneer  at  Station  6.  The  orange  coloration  was  replaced  by  brown  at 
Stations  7 through  11.  It  was  a brownish-orange  at  Station  12,  grading  to 
medium  brown  below.  Some  gravel  remains  of  an  old  dam  were  found  at  Station  9. 
Small  amounts  of  shell  were  observed  at  Stations  6 through  11. 

Detritus  was  present  at  all  stations,  but  was  abundant  only  at  Station 
6.  Oil  was  also  found  only  at  Station  6 and  formed  a film  on  the  water 
surface  when  sieving  macroinvertebrate  samples.  A mat  of  filamentous  algae 
covered  the  bottom  at  Stations  3 through  6 where  it  was  relatively  dense. 
However,  it  was  very  thin  at  Station  6.  A slender  spikerush  (Eleocharis 
acicularis?)  was  found  in  small  quantities  on  the  benthos  at  Station  IT. 

Dead  fish  were  present  in  the  sediment  at  Stations  1 through  5.  They 
were  partially  decomposed  and  almost  unidentifiable  but  were  probably  gizzard 
shad  and  appeared  most  abundant  at  Stations  1 through  3.  Burrows  of  the  may- 
fly larvae  Hexagenia  sp.  were  abundant  at  Stations  7 through  9 while  only  a 
few  were  found  at  Stations  10  and  11. 
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An  effort  was  made  to  sample  benthic  macroinvertebrates  at  comparable 
depths.  This  was  possible  at  all  stations  except  8,  9,  10  and  12.  Stations  8, 
9,  and  10  were  3-4  meters  deep,  and  Station  12  was  about  0.75  meters  deep.  The 
remaining  stations  were  1.3-2  meters  deep. 
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Description 


PHENOGR/VI  OF  PERIPHYTON,  ARTIFICIAL  SUBSTRATE, 
CULLED  1%,  MUTUAL  ABSENCE  IMPORTANT,  COPHENETIC 
CORRELATION  COEFFICIENT  .871. 
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COMPUTATIONAL  METHODS* 


Coimunitv 


Introduction 


Biotic  components  of  water  quality  are  generally  quantified  by  one- 
dimensional diversity  indices  when  single  samples  or  stations  are  examined, 
or  two-dimensional  coefficients  of  biotic  similarity  when  sample/sample, 
station/station,  or  species/species  comparisons  are  undertaken. 

Diversity  indices  are  mathematical  expressions  that  describe  the  distri- 
bution of  individuals  within  the  community.  There  are  a number  of  diversity 
expressions  in  use.  In  general,  maximum  diversity  exists  if  each  individual 
belongs  to  a different  species.  An  environmental  parameter  that  influences 
community  structure  will  also  modify  the  diversity  index.  In  cases  where 
environmental  stress  may  occur  (such  as  competition  among  species,  physio- 
chemical  limiting  factors,  or  pollution),  the  community  is  reduced  in  the 
number  of  species  present.  Frequently,  this  reduction  in  the  number  of 
species  is  accompanied  by  an  increase  in  the  number  of  individuals  of  the 
remaining  species,  especially  in  the  case  of  organic  pollution.  Environmental 
stress,  therefore,  tends  to  reduce  the  magnitude  of  diversity  indices.  One- 
dimensional diversity  indices  include  the  Shannon-Weaver  Species  Diversity, 
Evenness,  and  Simpson's  Index  of  Dominance. 

Coefficients  of  biotic  similarity  quantify  the  taxonomic  overlap  between 
two  samples  or  stations.  Most  of  these  coefficients  assume  values  between 
0 and  1,  where  a value  of  0 indicates  no  species  overlap,  and  a value  of  1 
implies  identical  species  composition.  Morisita's  Index  of  Faunal  Affinity 
and  the  Pinkham-Pearson  Index  of  Biotic  Similarity  are  measures  of  biotic 
similarity. 

In  this  study,  data  processing  subsequent  to  manual  taxonomic  identi- 
fication/confirmation was  executed  through  the  IBM  370/OS  system  at  the 
Northeast  Regional  Data  Center  of  the  State  University  System  of  Florida 
(NERDC).  Diversity  indices  and  coefficients  of  similarity  were  calculated 
by  proprietary  FORTRAN  IV  routines.  The  phenograms  were  generated  through 
application  of  the  NT-SYS  Numerical  Taxonomy  System  developed  by  Rohlf, 
Kishpaugh  and  Kirk  at  Stony  Brook  (1974). 


Shannon-Weaver  Species  Diversity  Index  (H 


The  Shannon  Weaver  Species  Diversity  Index,  H (Odum,  1971)  is  defined 
as : 

t n.  n. 
fr  N~ 


♦Procedures  outlined  and  examples  cited  were  part  of  a study  conducted  for  the 
U.S.  Army  Medical  Research  and  Development  Command. 
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where  = total  number  of  organisms  present  as  species  i^ 
t 

N = s n.  = total  number  of  organisms  present  in  the  sample 
i-1  ^ 

t - number  of  taxa  present  in  the  sample. 

H ranges  from  a minimum  of  0.0,  occurring  when  all  organisms  belong  to  the 
same  taxon  (no  diversity),  to  a maximum  of  In  N,  occurring  where  each 
organism  present  belongs  to  a unique  taxon  (maximum  diversity). 

The  Shannon-Weaver  Index  is  commonly  expressed  to  other  logarithmic 
bases,  especially  base  2 and  base  10,  and  is  easily  converted  by  the 
following  expression; 


H.  = 
basOj^  In  x 

Pinkham-Pearson  Index  of  Biotic  Similarity  (B) 

The  previously  discussed  index  (H)  quantifies  community  structure  with 
a sacrifice  of  taxonomic  integrity  important  to  paired  comparisons  between 
samples  or  stations.  Such  an  index  is  incapable  of  distinguishing  samples 
of  similar  gross  community  structure,  but  unlike  taxonomic  composition. 

That  is,  in  computation,  the  ^ th  species  of  one  sample  is  not  necessarily 
the  same  ]_  th  species  of  another  sample. 

This  insensitivity  to  taxonomic  overlap  is  surmounted  by  the  Pinkham- 
Pearson  Index  of  Biotic  Similarity,  B (Pinkham  and  Pearson,  1974)  defined 
as: 


R = i J Min  (”iA.  '^iB) 
t Max  (n!”^  n.g) 

where  t = number  of  taxa  considered 

n.^  = number  of  organisms  of  species  present  at  Station  A 
n^g  = number  of  organisms  of  species  j_  present  at  Station  B 
Min  (n^^,  n^.g)  = the  minimum  value  of  the  pair:  n.^,  n.g 
Max  (n^.^,  n^.g)  = the  maximum  value  of  the  pair:  n.^,  n^g 

Biotic  similarity  is  defined  only  for  a paired  comparison  between  two  samples 
or  stations.  If  two  samples  are  characterized  by  identical  taxonomic  overlap 
Ull  species  occur  in  identical  abundance),  the  calculated  index  assumes  a 
value  of  1.0  (maximum  similarity).  Two  samples  possessing  no  species  in  conmon 
share  an  index  of  0.0  (minimum  or  no  similarity).  The  number  of  species 
considered,  t,  may  include  only  those  species  observed  in  either  or  both  of 


the  two  samples,  or,  if  mutual  absence  is  deemed  important,  may  include 
species  not  necessarily  present  in  either  sample.  If  mutual  absence  is 
considered  important,  Min  (0,0)  s 1 and  Max  (0,0)  = 1 in  the  computation 
of  biotic  similarity. 

A biotic  similarity  index,  B',  between  species  may  be  defined  on 
spatial  and  numerical  occurrence  by  transposition  of  the  axes  in  the  pre- 
ceding expression  of  station  similarity; 


where: 


k = number  of  samples  or  stations  considered 

n.,  = number  of  organisms  of  species  1_  at  Station  i 

sJ  * 

nj2  = number  of  organisms  of  species  ^ at  Station  i 
Min  (nji*  11^2)  = the  minimum  value  of  the  pair:  nj-j. 
Max  (iiji,  nj2)  = the  maximum  value  of  the  pair:  nj^. 


"a2 

"j2 


This  index  likewise  ranges  from  0.0  (minimum  similarity)  to  1.0  (maximum 
similarity).  B'  may  possess  utility  for  grouping  species  according  to 
environmental  preference  or  pollution  tolerance  — that  is,  it  may 
delineate  "indicator  organisms". 


Phonograms 

The  quantification  of  similarity  between  paired  stations,  samples,  or 
species  by  any  of  the  previously-defined  coefficients  of  similarity  generates 
a diagonal  matrix  containing  PC  unique  elements,  where  PC  is  calculated  from 
the  expression  (Pinkham  and  Pearson,  1974): 

PC  . 

where:  PC  = number  of  unique  paired  comparisons 

S = number  of  stations,  samples,  or  species  being  compared. 

For  a study  comprising  only  25  stations,  a similarity  matrix  of  300  unique 
elements  is  produced.  Evaluation  and  presentation  of  such  a voluminous 
matrix  is  impractical  without  computer-aided  analysis  and  graphic  models. 

Algorithms  for  clustering  similarity  matrices  into  two-dimensional, 
hierarchic  relationships  have  been  developed  by  numerical  taxonomists  (Sokal 
and  Sneath,  1963).  A technique  frequently  invoked  by  ecologists  and  generally 
regarded  as  introducing  the  least  distortion  into  similarity  relationships 
is  the  sequential,  agglomerative,  hierarchic,  nonoverlapping  clustering 
method  (SAHN)  using  unweighted  pair-groups  with  arithmetic  averaging  (UPRMA), 
described  by  Sokal  and  Sneath  (1963).  The  product  of  this  procedure  is  a 
branched  diagram  termed  a phenogram  (or  dendrogram),  illustrated  in  Figure  R-1. 
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This  phenogram  was  generated  by  computer  using  the  NT-STS  Numerical  Taxonomy 
Package  (Rohlf,  Kishpaugh,  and  Kirk,  1974). 

The  horizontal  scale  or  abscissa  of  the  phenogram  is  graduated  in  the 
units  of  the  similarity  measure  upon  which  the  clustering  was  based  --  in 
this  case,  the  Pinkham-Pearson  Biotic  Similarity  Index.  Points  of  furcation 
(branching)  between  the  horizontal  stems,  representing  stations  or  groups  of 
stations  imply  that  the  similarity  between  the  two  streams  is  at  the  coefficient 
value  shown  above  the  branch  on  the  abscissa. 

The  magnitude  of  similarity  between  any  two  stations  represented  on  the 
phenogram  will,  in  general,  differ  from  the  corresponding  magnitude  given 
in  the  original  similarity  matrix.  This  arises  as  a consequence  of  the 
averaging  necessary  to  recursively  agglomerate  the  separate  stations  into  a 
single,  structured  set  containing  all  the  stations. 

The  degree  of  distortion  resulting  from  the  cluster  analysis  may  be 
quantified  by  the  cophenetic  correlation  coefficient,  r^pph,  defined  as  the 
product  moment  correlation  coefficient  computed  between  the  elements  of 
the  original  similarity  matrix  and  the  corresponding  indices  implied  by  the 
phenogram  (Sokal  and  Sneath,  1963).  High  values  of  '^(;oph(’^coph 
fewer  than  10  stations)  indicate  that  the  distortion  introduced  by  the 
clustering  procedure  and  depicted  by  the  phenogram  has  not  significantly 
masked  the  informational  content  of  the  original  similarity  matrix. 
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